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Executive  Summary 


Presented  in  this  report  is  detailed  information  describing  the 
need  for  a technical  base  for  air  conditioning  criteria.  The  report 
includes  results  of  a feasibility  study  for  evaluating  indoor  habitability 
conditions  for  air  conditioned  as  well  as  for  non-air  conditioned  build- 
ings during  summer  months.  In  the  feasibility  study  the  indoor  habita- 
bility condition  is  determined  by  calculation  using  a computerized  simu- 
lation technique  which  follows  a detailed  heat  transfer  analysis  of 
building-human  body  systems.  A concept  of  "Predicted  Indoor  Habitability 
Index"  (PIHI)  is  introduced  to  designate  an  index  for  describing  the 
indoor  thermal  environment.  PIHI  is  designed  as  an  index  that  relates 
the  calculated  indoor  temperature,  humidity,  air  velocity,  and  mean 
radiant  temperature  to  the  physiological  as  well  as  psychological  re- 
sponse of  occupants  in  such  a manner  that  the  short  term  as  well  as  long 
term  exposure  in  the  non-air  conditioned  buildings  can  be  evaluated. 

There  are  several  published  types  of  physiological  indices  (both  comfort 
and  discomfort).  A study  of  these  indices  with  respect  to  their  ap- 
plicability to  the  PIHI  concept  is  given. 
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In  the  introduction  the  need  for  an  air  conditioning  criteria  is 
given  together  with  the  manner  in  which  it  can  be  established.  This 
involves  a detailed  survey  of  building  data  and  weather  data  as  well  as 
a precise  thermal  simulation  of  buildings.  The  second  part  of  the  report 
describes  the  basic  mathematical  relationships  used  to  predict  indoor 
temperature  and  humidity  as  a response  of  the  building  to  climatic  factors. 
The  third  part  presents  progress  to  date  on  compiling  and  analyzing  build- 
ing data.  This  effort  is  essential  in  determining  air  conditioning  criteria 
and  entails  the  development  of  heat  transfer  characteristics,  air  leakage 
characteristics,  and  thermal  mass  characteristics  which  play  very  important 
roles  in  determining  the  indoor  conditions  of  any  building  unit.  An 
analysis  is  done  to  find  the  range  and  distribution  of  these  basic  parame- 
ters so  that  a meaningful  building  classification  useful  for  this  approach 
can  be  made.  The  third  part  will  be  the  subject  of  a separate  NBS  numbered 
report . 

In  Part  4,  a review  is  given  of  various  physiological  indices  currently 
available.  It  has  been  possible  to  include  very  recent  information  in 
this  part.  Part  5 gives  the  details  of  a comprehensive  feasibility  study 
of  a thermal  simulation  made  on  an  apartment  building  located  both  in 
Jersey  City,  New  Jersey  and  Macon,  Georgia.  Algorithms  for  many  of  the 
physiological  indices  were  used  in  conjunction  with  the  thermal  simulation 
to  give  hour  by  hour  values  of  these  indices  in  the  apartment.  An  examina- 
tion of  the  results  of  this  study  reveals  the  extent  and  duration  of  un- 
desirable indoor  conditions  when  this  apartment  is  not  air  conditioned. 

The  magnitude  of  computational  effort  appears  formidable  if  detailed 


ii 


simulation  calculations  are  to  be  carried  out  for  the  many  combinations 
of  building  parameters  and  climate  zones  in  the  United  States.  Therefore, 
an  attempt  was  made  to  obtain  statistical  correlations  between  the  indoor 
and  outdoor  conditions  for  this  apartment.  The  attempt  was  reasonably 
successful  and  a short-cut  technique  for  establishing  the  criteria  is 
indicated. 
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1.  Introduction 


At  the  present  time  in  the  United  States  government  agencies,  the 
decision  to  air  condition  a single  family  or  multifamily  dwelling  is 
usually  based  on  cost  and  general  location  of  the  housing  unit  (i.e., 
northern  or  southern  part  of  the  country).  The  United  States  Air  Force, 
for  example,  specifies  in  their  Manual  88-8  that  family  housing  located 
in  a region  where  the  wet -bulb  temperature  is  above  67  °F  for  1000  hours 
or  more  during  the  hottest  six  consecutive  months  of  the  year  should  have 
air  conditioning.  The  General  Services  Administration  uses  the  degree- 
day  concept  in  specifying  the  guidelines  for  air  conditioning  their  cars. 

If  there  are  more  than  700  cooling-degree-days  based  on  65  °F  in  a year 
(at  that  particular  location),  then  air  conditioning  should  be  installed. 

In  addition,  the  Office  of  Management  and  Budget  gives  a criteria  similar 
to  that  of  the  Air  Force  in  relation  to  the  air  conditioning  of  federal 
buildings.  For  lack  of  additional  information  on  family  housing,  the  Air 
Force  criteria  was  adopted  by  the  Department  of  Housing  and  Urban  Development 
and  is  present  in  their  four  volumes  "Guide  Criteria  for  the  Design  and  Eval- 
uation of  Operation  Breakthrough".  Using  this  as  a basis  for  the  decision 
on  air  conditioning  is  a substantial  improvement  over  no  criteria  at  all; 
however,  it  is  felt  that  other  factors  besides  weather  should  be  con- 
sidered. 
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One  of  the  most  important  items  that  govern  the  indoor  environment 
is  the  "thermal  performance"  of  the  enclosing  structure.  J.  F.  van  Straaten 
of  the  Building  Research  Institute  of  South  Africa  discusses  in  his  recent 
book  Thermal  Performance  of  Buildings  (1)  the  way  different  structures 
respond  when  subjected  to  periodically  varying  solar  conditions  and  en- 
vironmental temperatures.  The  overall  heat  transfer  coefficient  for  a 
wall,  which  gives  a direct  indication  of  the  resistance  to  heat  flow  under 
steady-state  conditions,  is  not  the  sole  indicator  of  the  heat  flow  to 
the  inside  environment  when  unsteady-flow  conditions  exist.  The  mass  of 
the  wall,  floor  and  ceiling,  which  governs  the  energy  storage  capacity, 
is  also  a pertinent  parameter  to  be  considered. 

Van  Straaten  reports  the  results  of  an  experimental  study  where  small 

2 

4 room  single  storey  structures  of  about  500  ft  floor  area  were  built, 
not  air  conditioned,  unoccupied,  and  instrumented  to  observe  the  indoor 
environment  as  a function  of  time.  Under  warm  arid  conditions  the  massive 
brick  buildings  proved  to  be  superior  in  minimizing  the  daily  variation 
of  the  indoor  temperature.  The  reason  is  that  the  "time  lag"  associated 
with  the  massive  structure  was  such  that  the  indoor  environment  remained 
sufficiently  cool  during  the  hot  part  of  the  day  (79  ®F  compared  with  84  °F 
outdoor  temperature  at  noon)  and  also  remained  sufficiently  warm  during 
the  night  (74  °F  compared  with  61  °F  outdoor  temperature  at  5:00  a.m.). 

The  light  timber  structures  with  essentially  no  "time  lag"  allowed  the 
indoor  temperature  to  increase  too  much  during  the  day  (93  ®F  compared 
with  87  ®F  outdoor  temperature  at  2:00  p.m.)  and  to  become  too  cool  during 
the  night  (62  °F  compared  with  61  °F  outdoor  temperature  at  5:00  a.m.). 
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On  the  other  hand,  van  Straaten  concludes  from  experiments  conducted 
in  Durban,  South  Africa  where  the  climate  is  warm  and  humid  and  the  daily 
outdoor  temperature  variation  is  considerably  less  even  though  the  daily 
average  temperature  is  approximately  the  same  and  solar  radiation  intensi- 
ties are  high,  that  massive  buildings  were  a distinct  disadvantage.  The 
walls  could  not  cool  down  sufficiently  during  the  night  to  provide  reason- 
ably comfortable  temperatures  in  the  enclosures  during  the  day.  Other 
interesting  results  of  van  Straaten *s  research  include  the  improvement 
in  performance  of  the  lightweight  structures  in  warm  and  dry  climates 
by  the  addition  of  concrete  floors  and  by  painting  the  exterior  of  the 
structure  white.  It  should  be  apparent  from  this  type  of  investigation 
that  climate  as  well  as  building  construction  plays  a significant  role 
in  determining  a need  for  air  conditioning. 

Obviously,  another  factor  to  be  considered  in  establishing  the  air 
conditioning  criteria  is  the  range  of  indoor  thermal  conditions  under 
which  people  remain  comfortable  or  not  overly  uncomfortable.  The  American 
Society  of  Heating,  Refrigerating  and  Air-Conditioning  Engineers  Standard 
55-66  is  very  specific  in  its  stipulation  of  thermal  comfort  conditions. 

In  a region  between  3 inches  above  the  floor  and  72  inches  above  the 
floor  and  at  least  2 feet  from  any  wall,  the  following  should  be  met  at 
all  times. 

1.  The  dry  bulb  temperature  shall  be  between  73  and  77  °F. 

2.  The  relative  humidity  shall  not  exceed  60  percent  and 
be  greater  than  20  percent. 

3.  The  air  motion  shall  not  exceed  45  fpm  and  be  greater 
than  10  fpm. 
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In  other  words,  if  the  indoor  thermal  environment  of  non-air  conditioned 
buildings  falls  into  this  condition,  it  is  obvious  that  the  building  needs 
no  air  conditioning  even  if  it  is  located  in  the  area  where  the  climatic 
criteria  indicates  the  need  for  it.  Because  of  the  inflexibility  of  the 
statement  "at  all  times" , it  may  not  be  appropriate  to  incorporate  this 
criteria  into  the  air  conditioning  design  criteria.  A better  evaluation 
of  the  tolerance  of  people  to  varying  indoor  conditions  above  and  below 
the  recommended  steady  state  comfort  conditions,  magnitude  and  frequency 
is  needed. 

There  are  also  other  basic  considerations  that  could  be  used  in 
establishing  an  air  conditioning  criteria.  For  example,  one  might  be 
able  to  maintain  reasonable  conditions  in  a space  by  a large  amount  of 
forced  or  natural  ventilation  yet  it  is  possible  that  it  would  not  be 
feasible  due  to  air  and/or  noise  pollution  problems.  It  might  be  cheaper 
to  design  completely  closed-in  structures  with  air  conditioning  than  more 
standard  type  buildings  without  air  conditioning.  Considerations  such  as 
these  will  not  be  dealt  with  initially  in  this  study. 

The  manner  in  which  the  criteria,  based  on  thermal  performance  alone, 
could  be  established  is  depicted  in  Figure  1.  All  combinations  of  selected 
weather  data  and  selected  building  data  are  combined  in  the  thermal  com- 
puter simulation,  the  output  of  which  is  the  space  conditions.  The  fre- 
quency and  duration  of  the  various  indoor  parameters  are  then  compared 
with  the  PIHI  to  determine  whether  the  space  should  be  conditioned  or 
not.  The  studies  completed  to  date,  discussed  in  Part  5,  show  the 
feasibility  of  such  an  approach. 
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It  should  also  be  noted  that  the  concept  of  determining  a need  for 
air  conditioning  based  on  building  thermal  performance  is  not  new. 

Givoni  in  his  latest  book  (2)  introduces  a "Building  Bioclimatic  Chart". 
It  is  essentially  a psychrometric  chart  on  which  he  has  drawn  an  en- 
closed region  indicating  thermal  comfort  conditions.  From  experience 
and  study,  three  additional  regions  have  been  drawn:  one  to  indicate 

the  range  of  outdoor  conditions  where  choice  of  building  materials  and 
design  features  (large  thermal  resistance,  large  thermal  mass,  etc.) 
make  comfort  conditions  inside  attainable  without  air  conditioning; 
a second  region  to  indicate  the  outdoor  conditions  where  efficient 
ventilation  through  the  structure  could  result  in  comfort  conditions;  and 
a last  region  showing  the  range  of  outdoor  conditions  where  mechanical 
cooling  would  be  required. 

In  a recent  article  of  the  ASHRAE  Journal  (3),  Nevins  and  McNall  of 
Kansas  State  University  discuss  a proposed  method  of  classifying  air  con- 
ditioning systems.  They  suggest  that  the  environment  of  a given  occupied 
space  can  be  evaluated  in  terms  of  the  human  response  to  that  environment 
and  that  an  appropriate  classification  or  "grade"  can  be  assigned  to  the 
system  serving  that  space  based  upon  the  functional  performance  in  the 
space.  Their  ideas  are  so  much  in  line  with  the  kind  of  analysis  pre- 
sented in  this  report  that  their  concluding  two  paragraphs  of  "Suggested 
Future  Action"  are  included  here: 

"We  suggest  that  a committee  be  formed  to  examine  the  various 
methods  of  establishing  a "percent  dissatisfied"  or  other 
functional  rating  method  for  an  environment  which  seems  most 
appropriate  in  the  light  of  all  present  data.  Such  a rating 
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method  must  take  into  account  practical  means  of  measuring  the 
environment  a 1 var iab le  s . 


Any  functional  specifications  must  then  be  met  with  practical 
equipment o Other  committees  must  assess  the  functional  per- 
formance of  systems  of  various  types.  The  system  in  this  case 
includes  not  only  the  air-conditioning  equipment  and  its  controls 
but  also  the  building,  its  construction  details  and  the  weather 
conditions  in  which  the  system  operates.  This  difficult  problem 
must  be  met  by  field  experiments,  that  establish  for  various  air- 
conditioning  equipment  and  control  combinations,  which  of  the  con- 
struction details  (window  sizes,  single  or  double,  etc.  and  wall 
*'U*'  values)  and  what  types  of  climates  (temperature,  wind,  solar 
exposure,  etc.)  provide  what  level  of  thermal  environment  in  the 
space  to  be  compared  with  the  corresponding  functional  "percent 
dissatisfied"  or  other  criteria.  When  all  this  is  accomplished, 
each  air-conditioning  system  and  controls  combination  could  be 
tabulated  with  ranges  of  climate  and  construction  details  which 
provide  various  classifications  of  thermal  (and  economical 
performance." 

Even  though  they  are  discussing  a different  end  product,  the  all-encom- 
passing analysis  suggested  is  very  similar  to  that  presented  in  this  report. 
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2.  Calculation  of  Space  Conditions  by  Building 
Thermal  Simulation  on  a Digital  Computer 


Numerous  papers  have  been  published  and  calculation  manuals  pre- 
pared in  the  area  of  heating/cooling  load  determination.  Very  few  of 
these,  however,  have  dealt  with  the  indoor  temperature  estimates  for 
the  condition  where  air  conditioning  is  non-existent  or  inadequate. 

The  room  temperature  prediction  is,  on  the  other  hand,  very  actively 
studied  in  Europe,  Japan,  Israel,  and  India  where  mechanical  cooling 
is  still  very  expensive  for  the  majority  of  dwellings  and  consequently 
the  design  of  a house  to  make  the  indoor  environment  habitable  is  ex- 
tremely important. 

The  heat  transfer  calculations  for  room  temperature  prediction 
are  similar  to  the  heating /coo ling  load  calculations;  the  former  are 
somewhat  more  complex  than  the  latter  because  they  require  exact  heat 
balance  calculations  for  the  room  air,  surrounding  walls  and  infiltrating 
outdoor  air.  Since  the  major  importance  is  the  hour  by  hour  room  tem- 
perature profile,  it  is  necessary  to  account  for  the  transient  heat 
conduction  and  storage  of  building  mass  as  affected  by  the  climatic 
parameters . 
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In  order  to  determine  temperature  and  humidity  of  non-air  conditioned 
rooms  responding  to  randomly  fluctuating  outdoor  climatic  conditions,  a 
computer  program  to  simulate  hourly  performance  of  building  heat  gain  and 
heat  storage  has  been  developed.  The  program  consists  of  various  sub- 
routines for  calculating  heat  gains,  which  are  similar  to  those  recommended 
by  the  ASHRA.E  Task  Group  on  Energy  Requirements  (4).  The  detail  of  these 
subroutines  are,  therefore,  not  given  in  this  report.  One  major  extension 
of  the  program  beyond  the  recommended  ASHEIAE  TG  subroutine  is  the  routine 
called  RMTMP,  which  solves  for  room  air  temperature.  General  mathematical 
expressions  of  this  routine  are  described  as  follows:  A general  expres- 

sion for  the  heat  exchange  at  an  interior  surface  of  the  room  is 


QK 


i,t 


M 


A. 

1 


^ = HC 


. (T  - T.  ) + S HD,  . (T,  ^ - T.  ^)  + RS.  ^ + RL.  ^ + RE.  ^ 
1 ^ a,t  i,t  ki  k,t  i,t  i,t  i,t  i,t 


(1) 


for  i = 1,  2 , . • . M 


where  QK.  ^ = rate  of  heat  conducted  into  surface  i at  the  interior 
i,t 


surface  at  time  t 


A.  = cross  sectional  area  of  the  surface 
1 


HC.  = convective  heat  transfer  coefficient  at  interior  surface  i 
1 

HD.  . = radiation  heat  transfer  factor  between  interior  surfaces 
ki 


i and  interior  surface  k 


T = room  air  temperature  at  time  t 
a,t 

T^  ^ = uniform  temperature  of  surface  i at  time  t 
T,  = uniform  temperature  of  surface  k at  time  t 
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RS.  ^ 


RL.  ^ 
i,t 


RE.  ^ 
i,t 


rate  of  solar  heat  coming  through  window  and  absorbed 
by  surface  i at  time  t 

rate  of  heat  radiated  from  lights  and  absorbed  by 
surface  i at  time  t 

rate  of  heat  radiated  from  appliances  and  occupants 
and  absorbed  by  surface  i at  time  t 


The  total  number  of  interior  surfaces  to  be  considered  for  this 

general  expression  is  designated  by  M and  it  varies  depending  upon  the 

complexity  of  the  room  construction.  If,  for  example,  each  of  four 

walls  are  made  of  four  different  segments  of  different  construction  and 

if  each  segment  of  the  wall  contains  a window  and  a door,  M will  be  2 x 

4x4+2  =34  with  the  assumption  that  ceiling  and  floor  are  of  uniform 

structure.  The  calculation  of  the  heat  conduction  QK.  is  performed  by 

1 , t 

the  response  factor  technique  for  heavy  construction,  such  as  exterior 
walls,  roof  and  floor  and  by  the  use  of  overall  heat  transfer  coefficient 
U (steady  state  calculation)  for  doors  and  windows  (lightweight  construc- 
tion). 

The  value  of  QK.  depends  upon  the  exterior  surface  temperature 
L , t 

T^  ^ at  the  opposite  side  of  the  given  interior  surface  by  the  following 
expression 


QK  .00  oo 

— 2 X.  . T.  ^ - S Y.  . T ^ . 

i,j  o,t-j 


A. 


(2) 


J-o 


1=0 


where  X.  . and  Y.  . (i  = 0,  1,  . . . oo)  are  the  response  factors  which  are 
1,1  1,1 

to  be  determined  as  functions  of  wall  (roof)  construction.  The  method 
to  calculate  the  response  factors  is  given  in  reference  (5). 
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On  the  other  hand,  the  following  expression  describes  the  heat 


transfer  process  at  the  exterior  of  the  surface: 


where 


QK 


2^ 


A. 

1 


QK 


o,t 


HC 

o 


HD 

,je 


RS 


o,t 


o,t 

DB. 


M' 

(T^  ^ - DB^ ) + 2 
o o , t t ^ 

i=i 


HD  (T  ^ - T ^)  - RS  ^ (3) 

OjA  o,t  i,t  o,t 


rate  of  heat  conducted  outward  at  the  exterior  side 
at  time  t 

convection  heat  transfer  coefficient  at  the  exterior 
side,  which  is  a function  of  wind  velocity 
radiation  heat  transfer  factor  between  the  exterior 
side  and  other  surface  i,  which  sees  the  surface*. 

Sky  and  ground  may  be  treated  as  one  of  these  I 
surfaces 

rate  of  solar  radiation  absorbed  at  the  exterior 
side  at  time  t 

temperature  of  the  exterior  surface  at  time  t 
outdoor  air  dry -bulb  temperature  at  time  t 


Again  in  Equation  (3)  the  conducted  heat  rate  QK  is  calculated  by  the 

o , t 

response  factors  for  the  heavyweight  structures  such  as  walls  and  roofs 
by 


QK  a. 

o < t _ ^ 

A.  ^ 


Y.  . T. 


i,j  i,t-j 


J=o 


- S Z.  . T ^ . 
ij  o,t-j 

j=o 


(4) 


where  Y.  . and  Z.  . are  the  response  factors  determined  from  the  con- 
struction  data.  For  the  lightweight  structures  such  as  doors  and  windows 


QK_. 


A. 

1 


= 


QK 


. T ) = U (T  - DB  ) 
A.  i,t  o,t  a,t  t 

1 

10 


(5) 


where  UT  is  the  overall  thermal  conductance  between  the  interior  and  ex- 


terior surfaces  and  U is  the  overall  heat  transfer  coefficient  between 
the  indoor  and  outdoor  air. 

In  addition  to  the  above  M equations  (equation  (1))  for  each  of  M 
interior  surfaces,  the  following  overall  heat  balance  equation  for  the 
room  has  to  be  solved  for  determining  the  room  air  temperature  T 


M 

S 

i=l 


T ^ + RS  ^ + RL  ^ + RE  ^ 
a,t  a,t  a,t  a,t 


(DB^  - T ^)  + pCG  (T  ^-T  ^)  = 0 

t a,t  ^ V v,t  a,t 


(6) 


where  RS  = rate  of  solar  heat  coming  through  windows  and  convected 
a ,t 

into  room  air  at  time  t 


RL  = rate  of  heat  from  light  convected  into  room  air  at  time 
a,  t 

t 


RE 


a,t 


P 

C 


rate  of  heat  from  appliances  and  occupants  convected 

into  room  air  at  time  t 

air  density 

air  specific  heat 

mass  flow  rate  of  outdoor  air  infiltrating  into  the 
room  at  time  t 

mass  flow  rate  of  ventilation  air  at  time  t 
outdoor  air  temperature  at  time  t 
ventilation  air  temperature  at  time  t 
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With  the  assumption  that  the  moisture  condensation  and  absorption 
over  the  room  surfaces  and  structures  is  small  the  room  humidity  is 
calculated  by  a simple  energy  balance: 


REL 


t ^ ^ <«o.t  - «a,t>  ^ ^ («v,t  - «a.t>  = ° 


(7) 


where  “ rate  of  latent  heat  released  into  the  room  air  by 

appliances  and  occupants  at  time  t 

\ = latent  heat  of  evaporation  of  water 

W^  ^ = humidity  ratio  of  outdoor  air  at  time  t 

W = humidity  ratio  of  ventilation  air  at  time  t 

v,t 

W = humidity  ratio  of  indoor  air  at  time  t 
a,t 

For  the  heat  dissipated  from  occupants,  which  are  included  in  Equations 
(1)  and  (6),  a linear  approximation  may  be  used  such  that 


Sensible  heat  release  of  time  t 
in  Btu/hr 


= 10  (100  - T 


Latent  heat  release  of  time  t 
in  Btu/hr 


= 10  (T  ^ - 60) 
a,t 


(8) 


where  T is  in  °F.  These  equations  are  valid  for  a resting  person 
a ,t 

(metabolic  heat  production  = 400  Btu/hr)  subjected  to  an  air  temperature 
between  60  ®F  and  100  °F.  Below  60  °F  all  the  heat  loss  can  be  considered 
as  sensible  and  above  100  °F  all  the  loss  can  be  considered  as  latent. 
l^.-ivchrv>iui't  r i c subroutines  permit  the  calculation  of  indoor  wet -bulb  tem- 
perature or  relative  humidity  with  the  room  air  humidity  ratio  W 

a,  t 

solved  for  from  Equation  (7). 
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For  Equations  (1)  and  (6)  the  radiative  and  convective  portions  of 


the  heat  emitted  from  lights  and  equipment  are  assumed  to  be 

Equipment,  RE^  ^ = 50%  of  total 

Appliances,  ’ 

and  People  M 

S RE.  • A.  = 50%  of  total 
i=i  ^ 


Light  RL  = 50%  of  total 

a , t 

M 

2 RL,  * A.  = 507o  of  total 
i=i  ^ 

The  radiation  portion  of  this  heat  gain  and  100%  of  the  net  solar  heat 
gain  (adjusted  for  the  shading  coefficient  coming  through  the  window) 
are  assumed  to  be  uniformly  distributed  and  absorbed  by  the  room  in- 
terior surfaces.  At  present  there  is  no  provision  to  include  the  effect 
of  shadow  over  the  exterior  surfaces  of  the  building  and  window  for  ad- 
justing the  solar  heat  gain. 

Although  it  is  well  known  that  air  infiltration,  G is  a function 

L,t 

of  outdoor  wind  velocity  and  air  leakage  opening  characteristics  of  the 

room  exterior,  accurate  relations  for  as  a function  of  wind  velocity 

L,  t 

and  building  structure  are  not  available.  Consequently  it  is  assumed, 
for  the  time  being,  that  the  air  infiltration  can  be  simulated  by  simply 
supplying  as  input  data  a specified  number  of  air  changes  per  hour.  The 
computer  program  presently  has  the  capability  of  changing  this  on  an 
hourly  basis  in  order  to  simulate  for  example  increasing  the  air  infil- 
tration at  night  for  purposes  of  more  effective  cooling. 
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Although  this  particular  computer  simulation  program  for  the  indoor 
temperature  and  hximidity  has  never  been  verified  by  experimental  obser- 
vation, many  similar  programs  of  the  same  basic  concept  have  been  known 

•k  j 

to  predict  the  indoor  temperature  to  within  2 °F~  . 

2.1  Indoor  Air  Velocity 

In  order  to  evaluate  indoor  conditions  for  human  comfort  or  physiologi- 
cal stress,  it  is  necessary  that  indoor  air  velocity  be  estimated.  The 
indoor  air  velocity  is,  however,  believed  to  vary  depending  upon  the  air 
leakage  rate,  ventilation  rate,  thermal  gradient  of  the  room,  basic  de- 
sign of  the  room,  and  any  local  air  moving  device  such  as  portable  fan. 

According  to  Givoni,  the  indoor  air  motion  could  vary  between  5 and 
60%  of  outdoor  wind  velocity  depending  upon  the  nature  of  the  air  passage 
resistance  through  the  room  (2).  Since  the  information  available  for 
predicting  indoor  air  velocity  as  function  of  temperature  distribution 
and  wind  velocity  is  very  scanty,  it  is  assumed  in  the  present  calculation 
that  the  indoor  air  velocity  is  less  than  20  ft  per  minute.  This  assumption 
should  yield  an  upper-bound  value  of  comfort  index  or  physiological  index 
for  summer  conditions  and  is  therefore  safer  to  use  than  to  overestimate 
the  cooling  due  to  elevated  air  velocities. 

* / 

~ Concurrently  a separate  investigation  is  taking  place  at  NBS  to  provide 
verification  of  the  computer  programs  on  several  HUD  Operation  Break- 
through buildings. 
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In  this  part,  a sample  format  for  the  compilation  and  presentation 
of  pertinent  building  data  relevant  to  the  thermal  performance  of  buildings 
is  shown.  Data  are  taken  from  the  proposals  of  the  22  participating 
housing  system  producers  (HSP)  in  the  HUD  Operation  Breakthrough  project. 
Since  the  22  HSP*s  involved  in  that  project  represent  the  major  housing 
system  producers  in  the  U.S.  building  industry,  it  is  felt  that  data  con- 
tained in  their  proposals  will  provide  adequate  and  meaningful  bases  for 
the  practice  and  the  state  of  the  art  in  housing  design. 

Up  to  the  present  time,  data  for  one-third  of  the  22  HSP*s  have 
been  compiled.  Because  of  the  different  and  numerous  ways  of  presenting 
their  design  by  the  contractors,  a computer  program  has  been  written  to 
arrange,  digest,  and  present  the  data  in  a unified  manner  for  easy 
reference.  Pertinent  parameters  associated  with  the  structural  components 
such  as  the  U-value,  weight  per  unit  area,  thermal  mass,  and  the  thermal 
time  constant  are  also  computed.  They  will  be  employed  for  the  purpose 
of  classification  of  the  structural  components  and  presented  together 
with  other  data  such  as  glass  area,  crack  length,  floor  area,  etc.  The 
physical  meaning  and  usefulness  of  these  parameters  will  be  discussed 
later.  Additional  parameters  will  be  computed  and  listed  if  and  when 
they  appear  to  be  of  significance  in  providing  information  in  the  clas- 
sification process.  A separate  report  will  be  prepared  giving  complete 
detail  of  this  process. 
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2 ,2  Mean  Radiant  Temperature 


Since  the  room  temperature  calculation  routine  is  providing  tempera- 
tures for  all  the  interior  surfaces,  it  is  an  easy  matter  to  determine 
the  mean  radiant  temperature  (MRT)  of  the  space.  Presently,  however, 
the  mean  radiant  temperature  routine  has  not  been  included  in  the  program 
and  it  is  assumed  that  MRT  is  the  same  as  the  air  temperature. 

3.  Classification  of  Building  Data 

One  of  the  most  important  things  that  govern  the  indoor  environment 
of  a building  is  the  "thermal  performance"  of  the  entire  structure.  Be- 
cause of  the  periodically  varying  solar  radiation  and  the  large  fluctuation 
in  external  environment  temperatures,  steady-state  heat  transfer  through 
the  enclosing  walls  into  the  interior  space  of  a building  rarely  occurs, 
and  the  heat  storage  effects  of  the  walls  and  the  internal  objects  must 
be  taken  into  account  in  order  to  estimate  accurately  the  cooling  load  or 
the  variations  of  the  indoor  environment.  The  overall  heat  transfer  co- 
efficient (U-value)  of  a wall,  which  gives  a direct  indication  of  the 
resistance  to  heat  flow  under  steady-state  conditions,  is  no  longer  the 
only  major  indicator  of  the  heat  flow  into  the  inside  environment  under 
transient  conditions.  The  mass  of  the  wall,  which  governs  the  energy 
storage  capacity,  is  also  a pertinent  parameter  to  be  considered. 
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3.1  Type  of  Building  Data  Compiled 


From  the  drawings  and  specifications  submitted  by  each  of  the 
housing  system  producers  in  the  HUD  Operation  Breakthrough  project,  the 
following  data  are  collected  for  each  of  the  different  types  and  sizes 
of  dwelling  units  contained  in  the  drawings: 

1,  Method  of  construction  of  the  exterior  wall,  party  wall, 
interior  partition,  roof,  ceiling/f loor , windows,  and  doors 

2,  Floor  area  for  each  of  the  rooms  in  a dwelling  unit 

3,  Type  of  walls  and  openings,  and  their  areas  for  each  of 
the  rooms  in  a unit 

4o  Crack  length  of  the  windows  and  doors 

5.  Interior  lighting  design 

6.  Natural  and  mechanical  ventilation  data 

7.  Shading  devices  for  the  window  areas 

3,2  Method  of  Building  Data  Presentation 

Data  collected  from  the  drawings  and  specifications  are  processed 
by  the  NBS  UNIVAC  1108  computer.  A sample  printout  for  a typical  dwelling 
unit  is  shown  on  the  enclosed  Tables  1 to  3 . It  should  be  pointed  out 
here  that  the  format  of  the  printouts  is  tentative,  and  may  be  revised 
to  a better  form  later  if  deemed  necessary. 


17 


Table  1 gives  a layer  by  layer  description  and  the  thermophysical 
properties  of  the  materials  used  in  the  construction  of  the  structural 
components.  The  order  of  the  layers  are  from  the  outside  of  a room 
toward  the  inside  for  the  side  walls  and  from  top  to  bottom  for  the 
roofs  and  floor /ceilings • The  thermal  property  values  are  taken  from 
the  HSP*s  drawings  and  specifications  from  the  1967  ASHRA.E  Guide  and 
Data  Book.  Very  often  data  for  the  specific  heat  of  certain 

materials  are  not  available,  and  a value  of  .200  Btu/lbm  °F  is  then 
assumed.  Since  the  value  of  does  not  vary  much  for  most  of  the 
building  materials  (of  the  order  of  .150  to  .400  with  a major  portion 
of  the  materials  in  the  lower  range  of  the  spectrum),  the  value  of  .200 
is  considered  a good  assumption. 

Table  2 gives  a room  by  room  tabulation  of  the  areas  of  the  floor, 
ceiling,  and  the  enclosing  walls  and  openings  of  a room  in  a dwelling 
unit.  The  four  sides  of  the  room  are  labelled  as  SI,  S2 , S3,  and  S4. 

The  orientation  of  the  unit  is  not  specified,  since  this  is  one  of  the 
factors  that  is  influenced  by  site  location,  terrain,  units  group  plan, 
etc.  as  much  as  by  thermal  performance  consideration.  For  a group  of 
attached  units  (single  family  attached,  multi-family  low  rise  or  high 
rise),  the  unit  at  the  ends  or  corner  of  the  group  is  chosen  since  it 
has  the  largest  amount  of  wall  area  exposed  to  the  outside  environment 
and  consequently  will  probably  result  in  the  most  severe  indoor  tempera- 
ture variation  or  the  largest  cooling  load  among  the  group.  Table  2 also 
gives  the  crack  lengths  of  the  windows  and  doors  which  are  required  in 
air  infiltration  specification.  Data  for  lighting  (as  internal  heat 
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source)  and  furniture  (as  internal  heat  sink)  are  not  available  at  the 
present  time,  and  will  be  incorporated  into  the  table  when  available. 

In  general,  data  contained  in  Table  2 are  required  by  the  computer  pro- 
gram developed  by  the  Environmental  Engineering  Section  of  the  NBS  Build- 
ing Research  Division  for  indoor  temperature  or  load  calculation  purposes 
It  should  be  pointed  out  that  to  simplify  the  data  read-off  procedure, 
certain  approximations  such  as  the  elimination  of  corners  and  the  neglect 
ing  of  closet  space  have  been  made.  The  data  in  Table  2 are,  therefore, 
not  an  exact  description  of  the  actual  floor  plan  in  the  HSP*s  drawings. 
However,  these  approximations  will  most  probably  not  affect  the  evalua- 
tion of  the  thermal  performance  of  the  dwelling  unit  in  any  significant 
way. 

Table  3 gives  the  values  of  various  relevant  parameters  associated 
with  and  affecting  the  thermal  performance  of  a building.  They  are  com- 
puted from  the  data  contained  in  Tables  1 and  2,  and  will  be  used  in  the 
classification  of  building  structural  components.  Briefly,  these  parame- 
ters are: 

1.  Floor  Area:  The  floor  area  represents  the  total  amount 

of  living  space  in  a dwelling  unit,  and  gives  an  indica- 
tion of  the  size  of  the  unit  among  the  same  type  of 
dwellings . 

2,  Ratio  of  Exterior  Wall  Area  to  Floor  Area:  This  ratio 

gives  an  indication  on  the  complexity  of  the  shape  of  the 
unit.  For  example,  a simple  rectangular  unit  will  have 

a lower  ratio  than  a more  complex  shaped  unit. 
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3,  Ratio  of  Window  Area  to  Exterior  Wall  Area:  This  ratio  in- 


4o 


5. 


dicates  the  relative  amount  of  glass  area  in  a unit  which 
has  a direct  bearing  on  the  amount  of  solar  radiation 
transmitted  directly  into  the  interior  space. 

U-value:  The  U-value  is  the  overall  heat  transfer  coef- 

ficient including  the  surface  air  film  resistance  (at  7.5 
mph  wind  outside  and  still  air  inside)  under  steady-state 
conditions.  This  is  the  major  factor  used  by  most  designers 
in  estimating  the  heating  and  cooling  load  of  a building. 

As  mentioned  previously,  under  unsteady-state  conditions 
the  U-value  is  no  longer  the  only  controlling  factor  in 
heat  flow  through  walls  since  it  does  not  take  into  account 
the  energy  storage  effect  of  the  walls.  However,  it  does 
give  an  indication  on  the  amount  of  resistance  the  wall 
offers  to  daily  average  heat  flow  and  therefore  acts  as 
one  of  the  parameters  in  determining  the  overall  heat 
transfer  through  the  walls. 

Thermal  Mass,  The  thermal  mass  is  a counterpart  of 

the  U-value  and  gives  an  indication  on  the  energy  storage 
capacity  of  the  walls.  It  is  computed  from  the  following 
equation: 


pC  L = 
P 


S 

i=l 


p.C  .L. 

pi  1 


(9) 


where  pi, 
thickness , 
structure. 


are  the  density,  specific  heat,  and 
respectively,  of  the  i^^  layer  in  a multilayer 
N is  the  total  number  of  layers  in  the  structure. 
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6.  Weight  per  unit  surface  area,  pL:  Because  of  the  small 


variations  in  the  magnitude  of  the  specific  heat  for  most 
building  materials,  the  heat  storage  capacity  of  a wall 
is  almost  directly  proportional  to  the  weight  per  unit 
surface  area  of  the  wall.  This  parameter  is  therefore 
closely  related  to  the  thermal  mass  of  the  wall.  It  has 
the  advantage  over  the  thermal  mass  in  that  a knowledge  of 
the  values  of  the  specific  heat  which  are  not  available 
for  certain  materials  is  not  required.  The  weight  per  unit 
area  of  a wall  is  computed  by 


7,  Thermal  Time  Constant  (TH,T.C,):  The  TH,T,C,  is  a function 

of  the  thermal  diffusivities  of  the  layers  in  a multilayer 
wall,  and  is  defined  as  the  heat  stored  per  unit  of  heat 
transmitted.  It  is  therefore  a combination  of  the  thermal 
mass  and  the  thermal  resistance  of  the  wall.  It  is  derived 
on  the  basis  of  the  analogy  between  the  heat  flow  represented 
by  a thermal  circuit,  and  the  time  constant  of  a R-C  electric 
circuit , and  is  given  by  (2 ) 


N 


(10) 


TH,T,C,  = S (Q/U)^ 


i 


(11) 


where 


+ L./2k.)  • (L.p.C  ) (12) 

1 1 11  pi 
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where  R is  the  resistance  of  the  surface  air  film  and 
os 

is  the  thermal  conductivity  of  the  material  of  the 
i-th  layer.  In  the  case  of  an  air  space,  sub- 

stituted by  the  resistance  of  the  air  space.  It  should 
be  pointed  out  that  the  TH.T.C.  given  by  the  above  equation 
also  takes  into  account  the  order  the  layers  of  materials 
are  arranged.  A change  in  the  order  of  the  arrangement  will 
change  the  rate  of  heat  flow  through  the  wall  under  transi- 
ent conditions.  This  is  because  that  the  outside  surface 
temperatures  of  the  walls  with  different  orders  of  arrange- 
ment in  their  layers  will  be  different  during  the  transient 
states.  For  example,  in  the  day  time,  the  surface  tempera- 
ture of  a wall  with  an  insulation  layer  next  to  the  external 
surface  layer  will  be  higher  than  a wall  constructed  of  the 
same  layers  of  materials  but  with  the  insulation  layer  placed 
further  away  from  the  external  surface  layer,  resulting  in 
more  heat  being  transferred  back  to  the  external  environment 
by  convection  and  radiation,  and  consequently  less  heat 
being  transferred  into  the  interior  space.  It  is  noted  that 
this  phenomena  is  not  accounted  for  by  either  the  U-value 
or  the  thermal  mass  where  the  order  of  arrangement  of  the 
layers  is  immaterial. 
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3.3  Planned  Work  in  the  Future  for  Building  Data  Processing 


After  the  completion  of  the  survey  of  building  data  on  the  22  housing 
system  producers,  a variety  of  selected  units  of  different  design  will  be 
used  as  input  to  the  computer  program,  and  the  resulting  indoor  tempera- 
ture variations  (not  air  conditioned)  or  cooling  load  (if  air  conditioned) 
under  the  same  external  environment  conditions  will  be  compared.  The 
influence  of  the  various  parameters  on  the  indoor  environment  will  be 
examined,  and  a classification  of  high  (or  large),  medium,  low  (or  small), 
or  various  ranges  of  numerical  values,  will  be  placed  on  those  parameters 
that  show  a definite  trend  in  relation  to  the  variations  in  the  indoor 
environment.  Additional  parameters  such  as  shading  factor,  ventilation 
rate  expressed  in  number  of  air  changes  per  hour,  and  color  of  the  ex- 
terior surfaces  will  also  be  examined.  A set  of  tables  based  on  the 
design  of  the  22  HSP’s  but  applicable  to  similar  types  of  design  by 
other  housing  producers  may  then  be  constructed.  The  number  of  the 
various  parameters  used  in  the  classification  will  be  kept  as  small  as 
possible  while  still  adequately  describing  the  structures  in  relation  to 
the  thermal  performance  of  a typical  dwelling  unit. 
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4o  Comfort  and  Physiological  Indices 


There  are  approximately  fifteen  prominent  indices  or  groups  of  ex- 
perimental data  that  have  been  presented  to  indicate  man’s  comfort  or 
discomfort  when  subjected  to  a certain  environment.  The  majority  are 
experimental  in  nature  (i.e,,  human  subjects  have  been  exposed  to  a 
specified  and  controlled  environment  and  asked  to  describe  their  feeling 
or  at  least  have  physiological  measurements  taken  on  them  during  the 
exposure).  However,  during  the  last  several  years  there  has  been  a 
great  deal  of  interest  in  describing  mathematically  the  thermal  inter- 
action between  man  and  his  environment,  thus  enabling  one  to  predict 
what  occurs  or  how  one  might  feel  in  given  surroundings.  It  is  generally 
accepted  that  some  form  of  skin  temperature  is  a good  measure  of  discom- 
fort in  cold  environments  whereas  sweat  rate,  skin  wettedness  or  some 
other  factor  related  to  sweat  production  is  the  best  indicator  of  dis- 
comfort in  warm  environments.  The  material  presented  here  will  con- 
centrate on  application  to  warm  environments.  The  remainder  of  this 
chapter  will  be  devoted  to  a discussion  of  the  most  widely  accepted 
of  these  indices . 
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4.1  Old  Effective  Temperature 


Undoubtedly  the  most  familiar  of  all  comfort  indices  is  the  American 

Society  of  Heating,  Refrigerating  and  Air-Conditioning  Engineers’  Effective 

Temperature.  The  scale  was  first  established  in  the  192 O’s  by  Houghten 

and  Yaglou  (6).  Subjects  wearing  clothing  having  an  unspecified  insula- 
*•  / 

tion  value”  were  passed  between  test  chambers  in  which  the  temperature, 
air  velocity,  and  humidity  were  maintained  at  specified  levels  and  were 
asked  to  describe  their  feeling  immediately  upon  entering  the  new  environ- 
ment. The  results  are  shown  in  Figure  2 where  lines  of  constant  thermal 
sensation  are  shown  plotted  on  a standard  psychrometric  chart.  The  ex- 
periments were  conducted  for  a large  range  of  air  velocities;  however, 

irk  / 

the  data  shown  here  is  for  an  air  velocity  of  less  than  20  fpm — . The 
lines  of  constant  thermal  sensation  were  arbitrarily  given  numerical 
values  corresponding  to  the  dry-bulb  temperature  at  which  the  line  inter- 
sected the  saturation  curve  (100%  relative  humidity).  If  the  dry-bulb 
temperature  were  a sole  indicator  of  thermal  comfort  according  to  this 
scale,  the  lines  would  run  vertically.  The  data  therefore  indicates  that 
* / 

” At  the  time  this  particular  index  was  developed,  the  analysis  of 

human  comfort  had  not  progressed  to  the  point  that  insulation  values 
for  clothing  ensembles  were  being  given. 

J 

— For  purposes  of  this  study,  the  data  for  all  the  indices  described 
will  be  for  very  nominal  air  velocities  (i.e.,  20  fpm  or  0,1  m/sec) 
since  higher  velocities  very  rarely  occur  in  residences. 
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as  the  relative  humidity  is  increased,  the  temperature  has  to  be  decreased 
to  maintain  the  same  feeling.  This  is  natural  since  cooling  by  sweat 
evaporation  is  restricted  somewhat  as  the  relative  humidity  increases. 

This  scale  originally  did  not  take  into  account  the  effect  of  radia- 
tion to  or  from  room  walls  where  their  temperature  is  maintained  at  values 
quite  different  from  the  room  air.  However,  a correction  was  proposed 
later  (7)  in  order  to  account  for  this  effect.  The  ET  scale  continues 
to  be  accepted  as  valid  in  the  very  hot  region  where  man’s  cooling  de- 
pends almost  entirely  on  his  sweat  evaporation  but  has  been  abandoned 
as  inaccurate  in  the  comfort  region  or  the  region  of  most  interest  to 
this  study. 

As  a result  of  the  apparent  disagreement  between  field  tests  and 
the  lines  of  constant  thermal  sensation  or  effective  temperature  at  the 
moderate  temperatures,  ASHRAE  undertook  an  additional  study  at  their  re- 
search laboratory  in  1960  and  the  results  are  shown  in  Figure  3.  In  this 
study  (8),  subjects  wore  light  indoor  clothing  and  were  seated  at  rest. 
Room  surfaces  were  held  at  the  same  temperature  as  the  room  air  in  all 
tests,  and  the  air  motion  in  the  space  was  20  fpm  or  less.  The  thermal 
sensation  votes  used  in  analyzing  the  test  data  were  those  which  were  cast 
after  approximately  three  hour  occupancy.  These  lines  show  the  conditions 
under  which  subjects  cast  votes  of  3,  4,  5,  or  6,  indicating  sensations 
of  slightly  cool,  comfortable,  slightly  warm,  and  warm  respectively.  It 
is  evident  from  Figures  2 and  3 that  the  old  effective  temperature  index 
predicts  a considerably  greater  effect  of  relative  humidity  on  the  thermal 
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sensation  feeling  than  do  the  lines  determined  by  Koch,  Jennings,  and 
Humphreys  (8).  In  defense  of  the  old  effective  temperature,  it  may  be  ex- 
plained that  the  two  groups  of  data  apply  to  different  conditions.  The 
effective  temperature  lines  indicate  a person *s  feeling  immediately  after 
entering  the  conditioned  space,  while  the  lines  of  Figure  3 are  for  people 
who  have  been  in  the  conditioned  space  for  approximately  three  hours  and 
are  in  equilibrium  with  the  environmental  conditions. 

4,2  Resultant  Temperature 

In  1948,  Missenard  (9)  proposed  an  index  called  the  resultant  tem- 
perature. Its  development  was  motivated  by  the  fact  that  the  effective 
temperature  was  not  based  on  experiments  where  thermal  equilibrium  was 
achieved  between  the  human  body  and  the  environment.  The  experiments  on 
which  the  resultant  temperature  was  based  were  conducted  with  about  six 
subjects  in  a psychrometric  chamber  with  a duration  of  exposure  greater 
than  was  used  in  the  experiments  for  the  effective  temperature.  This 
index  is  identical  in  structure  with  the  effective  temperature  and  is 
shown  in  Figure  4 for  clothed  subjects  exposed  to  an  air  velocity  of 
0.1  m/sec.  The  resultant  temperature  predicts  a much  less  effect  of 
thermal  sensation  on  humidity  than  did  the  effective  temperature.  How- 
ever, the  results  still  indicate  more  of  a dependence  on  humidity  in  the 
moderate  temperature  range  than  do  the  comfort  votes  of  Koch,  et.al. 
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4.3  Kansas  State  Index 


Perhaps  the  most  definitive  experimental  study  on  thermal  comfort 
is  one  conducted  on  over  1600  college  students  at  Kansas  State  University 
over  the  past  several  years  (10,  11).  The  students  were  exposed  (over 
a 3 hour  period)  to  a uniformly  heated  or  cooled  environment  (wall  tem- 
peratures identical  to  the  air  temperature)  with  moderate  air  velocities 
and  asked  to  vote  on  a scale  from  1 to  7 which  ranged  from  cold  to  hot: 

1.  cold 

2.  cool 

3.  slightly  cool 
4o  comfortable 

5.  slightly  warm 

6 . warm 

7 . hot 

All  subjects,  males  and  females  alike,  were  clothed  in  cotton  twill 

shirts  and  trousers.  The  shirts  were  worn  outside  of  the  trousers. 

The  subjects  wore  the  underwear  they  had  on  when  reporting  for  this 

test,  which  usually  consisted  of  cotton  undershorts  or  jockey  shorts 

for  the  men  and  brassieres  and  underpants  for  the  women.  All  other 

clothing,  including  T-shirts,  girdles,  slips,  etc.,  were  removed.  The 

subjects  wore  cotton  sweat  socks  without  shoes.  The  net  insulating 

value  for  the  clothing  ensemble,  as  computed  by  means  of  an  electrically 

^ / 

heated  copper  manikin  (12)  has  been  determined  to  be  0.6  do—  v/hich 
has  been  generally  accepted  as  the  standard  clothing  for  thermal  comfort 
studies . 


1 do  represents  a measured  resistance  to  heat  transfer  of  the  clothing 
ensemble  of  0.18  sq.  m hr  °C/Kcal  (0.88  sq.  ft  hr  °F/Btu  or  equivalent  of 
an  overall  heat  transfer  coefficient  = 1.13  Btu  per  (hr)  (sq.  ft)  (F)) 
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Even  though  the  results  are  available  for  men  and  women  separately 
and  also  for  the  votes  taken  at  each  half  an  hour  for  the  entire  3 hours 
of  exposure,  the  data  shown  in  Figure  5 is  an  average  for  the  votes  of 
men  and  women  taken  together  at  the  end  of  the  3 hours.  Rohles  (13) 
discusses  the  fact  that  a single  line  on  a temperature  grid  to  represent 
a "comfortable"  thermal  sensation  provides  little  in  the  way  of  practical 
information  unless  some  type  of  limits  surround  the  line  or  some  border 
is  presented  which  differentiates  the  "slightly  cool"  conditions  from 
the  "comfortable"  conditions  and  in  turn  the  "comfortable"  dimension 
from  the  limits  of  the  "slightly  warm"  conditions.  Consequently,  the 
data  was  reanalyzed  and  the  modal  comfort  envelopes  shown  in  Figures  6 
and  7 were  presented.  As  long  as  the  temperature  and  humidity  gave  a 
condition  within  the  70%  envelope,  at  least  70?o  of  the  subjects  voted 
4 (comfortable).  The  same  interpretation  goes  with  the  60%  envelope. 

It  should  be  pointed  out  that  over  907o  of  the  subjects  voted  either 
3,  4,  or  5 (slightly  cool,  comfortable,  or  slightly  warm)  when  the 
space  condition  was  anjA^/here  within  the  707o  envelope. 


29 


4,4  Predicted  Mean  Vote  and  Predicted  Percentage  of  Dissatisfied 


One  of  the  most  significant  attempts  to  predict  analytically  the 
thermal  interaction  between  man  and  his  environment  is  presented  by 
P.  0.  Fanger  of  the  Technical  University  of  Denmark  in  his  recent  book 
Thermal  Comfort  (14).  He  begins  by  pointing  out  the  most  important 
variables  which  influence  the  condition  of  thermal  comfort: 
lo  activity  level  (heat  production  in  the  body) 

2,  thermal  resistance  of  the  clothing  (clo-value) 

3.  air  temperature  (dry-bulb) 

•k  / 

4o  mean  radiant  temperature”  (to  account  for  the  radiation 
transfer  between  the  subject  and  the  surroundings) 

5.  relative  air  velocity 

6.  water  vapor  pressure  in  ambient  air 

The  basis  of  his  analysis  in  the  so-called  general  comfort  equation  which 
defines  all  combinations  of  the  variables  which  will  create  thermal  com- 
fort. Its  specific  form  is  as  follows: 


M - E 


diff 


- E - E 

sw 


res 


-D  = K = R + C 


(13) 


•k  / 

— Mean  radiant  temperature  is  defined  as  the  uniform  surface  temperature 
of  an  imaginary  black  enclosure  with  which  man  (also  assumed  a black 
body)  exchanges  the  same  heat  by  radiation  as  in  the  actual  environ- 
ment . 
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where 


M = the  internal  heat  production  in  the  human  body 
= the  energy  loss  by  water  vapor  diffusion 
through  the  skin 

E = the  energy  loss  by  evaporation  of  sweat  from 

the  surface  of  the  skin 

E = the  latent  respiration  heat  loss 
res 

D = the  dry  respiration  loss 

K = the  heat  transfer  from  the  skin  to  the  outer 

surface  of  the  clothed  body  (conduction  through 
the  clothing) 

R = the  heat  loss  by  radiation  from  the  outer  surface 
of  the  clothed  body 

C = the  heat  loss  by  convection  from  the  outer  surface 
of  the  clothed  body 

The  double  equation  expresses  that  the  internal  heat  production  M 
minus  the  heat  loss  by  evaporation  from  the  skin  ^sw^ 

respiration  (E  + D)  is  equal  to  the  heat  conducted  through  the 

1T0  S 

clothing  K and  dissipated  at  the  outer  surface  of  the  clothing  by 

radiation  and  convection  (R  + C).  It  is  assumed  that  the  evaporation 

corresponding  to  E and  E,  takes  place  at  (or  underneath)  the  skin 

sw  dirt 

surface. 
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After  formulating  expressions  for  each  component  in  the  above  equa- 
tion (see  Appendix  A)  and  substituting  them  into  Equation  13,  two  algebraic 
equations  result  in  the  following  9 unknowns: 

M - (defined  previously) 

Tg  - skin  surface  temperature 

- water  vapor  pressure  in  the  ambient  air 

E - (defined  previously) 
sw 

- ambient  air  dry  bulb  temperature 

" clothing  surface  temperature 

I - - insulation  value  of  the  clothing 
ci 

^MRT  " radiant  temperature  of  the  environment 

V - the  relative  air  velocity 

Six  of  the  unknowns  (M,  P^,  T^,  ^MRT^  considered  as 

controllable  or  specified  for  any  particular  subject --environment  situa- 
tion. Two  of  the  remaining  unknowns  are  eliminated  by  what  may  be  con- 
sidered Fanger*s  most  questionable  assumption.  During  studies  at  Kansas 
State  University  on  college  students,  data  was  obtained  for  mean  skin 
temperature  and  evaporative  heat  loss  of  the  subjects  while  they  reported 
being  comfortable  conducting  various  activities  (different  M*s).  E 

sw 

and  were  plotted  individually  as  a function  of  the  metabolic  production 
and  even  though  the  data  was  rather  widely  scattered,  linear  equations 

T = fct  (M)  (14) 

s 

and 

E = fct  (M)  (15) 

sw  ' 
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were  determined  by  regression  analysis.  The  procedure  then  used  was  to 

specify  5 of  the  6 controllable  or  specified  quantities  (M,  P^,  T^, 

V,  and  T ),  and  use  the  4 equations  (13,  14  and  15)  to  determine  the 
rlRl 

remaining  four  unknowns,  one  of  which  is  the  sixth  of  the  controllable 
parameters  that  would  give  thermal  comfort.  The  results  are  presented 
in  the  form  of  comfort  charts. 

Having  determined  the  conditions  for  optimum  thermal  comfort, 

Fanger  recognized  the  need  for  being  able  to  evaluate  an  existing  room 

climate  and  quantify  its  deviation  from  the  comfort  condition.  A similar 

scale  to  the  one  used  in  the  Kansas  State  studies  was  chosen; 

-3  cold 
-2  cool 

-1  slightly  cool 
0 neutral 
+1  slightly  warm 
+2  warm 
+3  hot 

The  numerical  values  however  are  lower  by  4.  A scale  is  thus  obtained 
which  is  easier  to  remember,  as  it  is  symmetrical  around  the  zero  point, 
so  that  a positive  value  corresponds  to  the  warm  side  and  a negative 
value  to  the  cold  side  of  neutral.  He  then  proposed  that  the  thermal 
sensation  values,  to  be  called  the  Predicted  Mean  Vote  (PMV),  would  be  a 
function  of  the  thermal  load  of  the  body,  defined  as  the  difference  be- 
tween the  internal  heat  production  and  the  heat  loss  to  the  actual  en- 
vironment for  a man  hypothetically  kept  at  the  comfort  values  of  the  mean 
skin  temperature  and  the  sweat  secretion  of  the  actual  activity  level. 

PMV  = fct  (M  - E,...  - E - E - D - R - C)  (16) 

diff  sw  res 
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The  nature  of  Equation  16  enabled  Fanger  to  use  the  same  component  ex- 
pressions as  in  the  general  comfort  equation.  The  particular  form  of 
the  function  was  determined  from  analyzing  the  Kansas  State  studies 
once  again  as  well  as  similar  experiments  conducted  at  the  Technical 
University  of  Denmark.  The  exact  expression  corresponding  to  Equation 
16  is  given  in  Appendix  A. 

Values  of  the  predicted  mean  vote  (PMV)  are  shown  on  the  psychrometric 

chart  in  Figure  8.  The  data  was  calculated  assuming  a resting  metabolic 

2 

rate  (M  = 50  kcal/hr  m ) , an  air  velocity  of  0.1  m/sec,  a mean  radiant 
temperature  equal  to  the  ambient  air  temperature,  and  a clothing  insula- 
tion value  of  0.6  do.  In  Figure  9 the  results  are  compared  with  the 
experimentally  determined  KSU  comfort  vote.  As  can  be  seen  the  agreement 
is  excellent  for  the  comfort  line  (KSU  = 0,  PMV  =4)  but  not  so  good  for 
environmental  conditions  outside  of  the  comfort  range. 

An  advantage  of  an  analytical  model  such  as  this  is  the  capability 
of  observing  the  effect  of  one  of  the  six  controllable  parameters 
without  having  to  conduct  extensive  experiments*  For  example,  it  may 
be  unrealistic  to  determine  an  air  conditioning  criteria  using  sub- 
jects that  wear  the  standard  0.6  do  clothing  ensemble.  Certainly 
in  hot  weather  people  would  be  more  inclined  to  wear  lighter  clothing 
such  as  shorts  and  an  open  neck  shirt  with  short  sleeves.  Data  for  this 
type  of  clothing  ensemble  (do  = 0.25)  is  shown  in  Figure  10  and  compared 
with  the  standard  Kansas  State  clothing  ensemble  in  Figure  11.  The  com- 
parison shows  what  one  might  intuitively  expect.  In  the  ligher  clothing, 
subjects  would  report  comfort  at  4 or  5 °F  warmer  temperature.  It  would 
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be  expected  that  they  would  vote  cool  at  a somewhat  higher  temperature 
also  (approximately  6 ®F).  However,  both  subjects  would  vote  hot  at 
approximately  the  same  temperature  and  relative  humidity. 

Instead  of  just  giving  the  predicted  mean  vote  as  an  expression  for 
the  thermal  environment,  Fanger  felt  it  might  be  even  more  meaningful  to 
state  what  percent  of  persons  can  be  expected  to  be  decidedly  dissatisfied 
since  this  can  readily  be  interpreted  by  both  the  engineer  and  the  layman. 
This  same  concept  has  been  proposed  most  recently  by  Nevins  and  McNall 
(3).  After  choosing  the  limits  of  -2  and  +2  for  being  dissatisfied  on 
the  cold  and  hot  side  respectively,  the  American  and  Danish  experiments 
already  mentioned  were  again  analyzed  to  determine  the  percent  of  the 
subjects  voting  cool  or  below  and  warm  or  above  at  each  environment  con- 
dition. The  same  relationship  was  assumed  to  hold  between  PMV  and  the 
new  factor,  predicted  percentage  of  dissatisfied  (PPD).  The  exact  rela- 
tion is  given  in  Appendix  A and  some  results  using  the  PPD  are  shown  in 
the  next  chapter. 

4.5  New  Effective  Temperature 

The  other  major  analytical  contribution  to  the  field  of  thermal 
comfort  to  be  mentioned  in  this  paper  is  a transient  heat  transfer 
model  presented  by  Gagge,  Stolwijk,  and  Nishi  (15).  The  model  differs 
most  significantly  from  Fanger *s  in  the  fact  that  physiological  responses 
have  been  included  such  as  blood  flow  restriction  and  dilation  to  and 
from  the  skin  and  regulatory  sweating  proportional  to  a deep  body  tem- 
perature deviation  from  some  predetermined  set  point.  The  human  body 
is  pictured  as  composed  of  two  distinct  parts;  a central  core  and  a 
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skin  shell.  An  energy  balance  is  written  for  each  part 


the  rate  of 

the  net  rate 

the  net  rate 

energy  in- 
crease of 
the  skin  shell 

of  energy 
transfer  from 
the  core  to 
the  skin 

+ 

of  energy 
transfer  from 
the  environ- 
ment to  the 
skin 

and 


the  rate 

the  rate 

the  net  rate 

— — 

the  net  rate  of 

of  energy 

of  energy 

of  energy 

energy  transfer 

increase 

generation 

+ 

transfer  from 

+ 

from  the  environ- 

of the 

(metabolic 

the  skin  to 

ment  to  the  core 

central 

production) 

the  core 

- 

core 

in  the  core 

or  in  equation  form 


dT 


^sk^sk 


sk 


d'^  ^con  ^bl  ^diff 


E - R - C 
rsw 


(16) 


and 


where 


dT 

M C -7^  = M - Q 
cr  cr  d'^  con 


E 

res 


M , = mass  of  the  skin  shell 
sk 

= specific  heat  of  the  skin  shell 
= temperature  of  the  skin  shell 


(17) 


M = mass  of  the  central  core 
cr 


C = specific  heat  of  the  central  core 
cr 

= temperature  of  the  central  core 
t = time 
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Q = heat  transferred  by  conduction  from  the  core  to 
con 

the  skin  (a  function  of  T and  T , ) 

cr  sk' 

Qbi  ~ energy  transferred  by  blood  flow  from  the  core 

to  the  skin  (a  function  of  T and  T , ) 

cr  sk 

= the  energy  loss  by  water  vapor  diffusion  through 
the  skin 

Ersw  “ energy  loss  by  evaporation  as  a result  of 

regulatory  sweat  production  (a  function  of 

T and  T , )l/ 
cr  sk 

R = the  heat  loss  by  radiation  from  the  skin  to  the 
environment 

C = the  heat  loss  by  convection  from  the  skin  to  the 
environment 

M = metabolic  heat  production  in  the  core 
E^^^  = the  latent  respiration  loss  from  the  lungs 

The  form  of  the  differential  equation  implies  the  lumped  capacity  analysis 
or  in  other  words,  the  temperature  of  the  skin  and  core  are  assumed  uni- 
form throughout  at  every  instant.  The  model  is  no  doubt  unrealistic  in 
this  respect,  since  the  temperature  can  and  does  often  vary  widely  within 
the  body. 

* / 

•“  This  particular  component  is  extremely  important  in  the  performance 
of  the  model  yet  the  expression  doesn't  agree  with  all  physiological 
observations  reported  in  the  literature.  A discussion  of  this  dis- 
crepancy is  given  in  Appendix  B. 

I 
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After  formulating  expressions  for  the  various  energy  components 

and  choosing  appropriate  numerical  values  for  the  various  parameters 

(see  Appendix  A),  the  two  differential  equations  are  solved  simultaneously 

(numerically  on  a digital  computer)  to  predict  the  time  variation  of  the 

core  and  skin  temperatures  when  the  controllable  parameters  mentioned 

previously  are  specified.  The  model  predicts  a response  that  levels 

off  at  the  end  of  approximately  one  hour.  Once  T and  T , have  been 

cr  s K 

determined,  it  is  possible  to  work  back  through  the  analysis  and  calcu- 
late the  regulatory  sweat  production  (at  the  end  of  one  hour),  associated 
with  the  given  environmental -human  conditions.  Mass  transfer  theory  is 
then  used  to  calculate  the  maximum  possible  evaporation  cooling  for 
given  skin  temperature  and  ambient  temperature  and  vapor  pressure. 

The  ratio  of  calculated  regulatory  sweat  evaporation  to  the  maximum 
possible  is  given  the  notation  of  skin  wettedness.  It  is  then  posulated 
(as  a result  of  experimental  observation)  that  in  warm  environments, 
various  conditions  giving  the  same  skin  wettedness  will  also  cause  the 
same  thermal  sensation  or  feeling  of  warmth  or  coolness.  This  particular 
point  is  debatable  and  is  the  subject  of  much  discussion  in  Appendix  B. 
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Lines  of  constant  skin  wettedness  at  the  end  of  one  hour’s  exposure 

2 

are  shown  in  Figure  12  for  a resting  subject  (58.2  watts/m  ) with  the 
standard  clothing  (0.6  do)  subjected  to  a small  air  velocity.  This 
model  has  been  proposed  to  ASHRAE  and  accepted  as  properly  predicting 
the  lines  of  constant  thermal  sensation.  Consequently  the  lines  in 
Figure  12  are  labelled  as  lines  of  New  Effective  Temperature.  Actually 
the  model  predicts  zero  regulatory  sweating  at  77  °F  and  507o  R.H.  or  in 
other  words,  the  77  ET  lines  corresponds  to  a skin  wettedness  of  0. 
Effective  temperature  lines  below  this  temperature  are  simply  drawn  as 
parallel  to  the  77  ET  line.  The  numerical  values  associated  with  the 
lines  were  determined  by  using  the  dry  bulb  temperature  at  the  intersection 
of  the  constant  skin  wettedness  line  and  the  507.  relative  humidity  curve. 
This  in  contrast  to  the  original  effective  temperature  where  experimentally 
determine  lines  of  constant  thermal  sensation  were  given  numerical  values 
corresponding  to  the  dry  bulb  temperature  where  the  lines  intersected 
the  saturation  or  100%  relative  humidity  curve.  The  labelling  of  the 
new  effective  temperature  seems  more  reasonable  since  the  values  will 
more  closely  correspond  to  dry  bulb  temperatures  experienced  in  everyday 
living.  However,  the  confusion  that  may  result  after  publication  of 
the  1972  volume  of  the  ASHRAE  Handbook  of  Fundamentals  can  be  seen  by 
observing  Figure  13.  The  numerical  values  of  the  two  scales  are  very 
widely  separated.  Since  the  old  effective  temperature  was  considered 
standard  for  almost  50  years,  many  guidelines  and  comfort  limits  will 
have  to  be  reevaluated  in  terms  of  the  new  scale. 
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By  observing  the  new  scale,  one  sees  that  the  dependence  of  thermal 
sensation  on  relative  humidity  at  the  higher  temperatures  is  approximately 
the  same  as  with  the  old  effective  temperature.  This  is  quite  acceptable 
since  the  old  effective  temperature  scale  was  considered  accurate  at 
these  temperatures.  However  at  temperatures  near  the  comfort  zone,  the 
new  model  predicts  a much  less  dependence  on  relative  humidity  which  is 
in  good  agreement  with  the  most  recent  studies.  These  two  factors  no 
doubt  contributed  significantly  to  the  new  scale's  acceptance  by  ASHRAE, 
Published  with  the  new  scale  in  the  1972  Handbook  will  be  a comfort  zone 
and  a broad  comfort  zone  shown  in  Figures  14  and  15  respectively. 

4.6  Heat  Stress  Index 

The  heat  stress  index  was  developed  by  Belding  and  Hatch  at  the 
University  of  Pittsburgh  and  first  published  in  1955  (16).  The  postula- 
tion that  forms  the  basis  of  their  index  is  as  follows:  the  metabolic 

production  of  the  body  minus  the  energy  that  can  be  given  off  by  radia- 
tion and  convection  to  the  surroundings  equals  a quantity  of  energy  that 
must  be  given  off  by  sweat  evaporation  in  order  to  maintain  thermal 
equilibrium  with  the  environment.  Consequently  the  expression  for  the 
required  evaporative  cooling  is 

E = M - R - C (18) 

req 

As  mentioned  previously,  for  a given  skin  temperature  and  ambient  air 
temperature  and  vapor  pressure,  mass  transfer  theory  can  be  used  to 
calculate  the  maximum  evaporative  cooling  possible.  The  heat  stress 
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index  is  defined  as  being  equal  to  the  ratio  of  the  two  quantities 

E /E  multiplied  by  100,  Belding  and  Hatch  assumed  a constant 
req  max 

skin  temperature  of  35  ®C  and  used  the  most  up  to  date  (at  that  time) 
expressions  for  radiation  and  convection  heat  loss  from  the  body  (see 
Appendix  A)  to  calculate  values  of  the  index  and  present  them  in  nomo- 
graphs as  functions  of  the  metabolic  production  and  environmental  con- 
ditions, In  addition,  physiological  and  hygienic  implications  of  8 hour 
exposures  to  various  heat  stresses  were  given  and  are  shown  in  Table  4. 

Calculations  have  been  made  using  the  original  expressions  for 
radiation  and  convection  losses  and  maximum  evaporative  cooling  from 
the  body  for  a resting  person  subjected  to  a uniform  environment  with 
an  air  velocity  of  20  fpm.  The  results  are  shown  in  Figure  16  in  the 
form  of  lines  of  constant  heat  stress  index.  Notice  that  using  these 
expressions  sweat  evaporation  to  maintain  body  equilibrium  is  not  re- 
quired until  the  environment  reaches  approximately  82  °F,  This  is  in 
definite  contrast  to  Gagge*s  effective  temperature  model  which  shows 
regulatory  sweat  production  starting  at  approximately  5 °F  lower.  The 
expressions  were  reevaluated  and  updated  in  1966  by  McKarns  and  Brief 
(17).  ASHRAE  plans  to  recommend  in  their  1972  Handbook  of  Fundamentals 
the  calculation  of  the  index  using  the  same  expressions  for  R,  C and 

E that  were  used  in  the  new  effective  temperature  model.  In  addition, 
max 

expressions  for  dry  and  latent  heat  losses  from  the  lung  passages  have 
been  subtracted  from  the  right  side  of  Equation  18  and  the  results  of 
new  calculations  are  shown  in  Figure  17  and  compared  with  the  original 
lines  of  constant  heat  stress  in  Figure  18.  The  latest  computations 
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seem  to  be  reasonable  since  the  evaporative  cooling  as  a result  of 
sweat  production  is  seen  to  be  necessary  at  a temperature  as  low  as 
74  °F.  Figure  19  shows  that  a heat  stress  index  of  10  which  corresponds 
to  mild  heat  strain  agrees  quite  well  with  the  upper  side  of  the  modal 
comfort  envelopes  as  proposed  by  Rohles  (13). 

As  the  new  computations  for  heat  stress  index  were  being  carried 
out , special  note  was  made  of  the  addition  of  expression  for  dry  and 
latent  heat  losses  from  the  lung  passages  to  Equation  18.  This  seems 
quite  reasonable  since  this  cooling  occurs  naturally  and  would  certainly 
assist  in  maintaining  body  thermal  equilibrium.  However , since  latent 
cooling  due  to  water  vapor  diffusion  through  the  skin  also  occurs 
naturally  the  computations  were  carried  out  once  again  subtracting  this 
factor  ) from  the  right  side  of  Equation  18.  The  results  are 

shown  in  Figure  20  and  it  can  be  seen  that  they  predict  evaporative 
cooling  due  to  regulatory  sweating  necessary  at  approximately  77  °F 
(50%  relative  humidity)  which  agrees  almost  identically  with  the  new 
effective  temperature  mode 1 • The  contribution  that  the  three  "natural" 
cooling  processes  make  can  be  seen  in  Figure  2 1 where  on  the  one  hand 
only  the  new  expression  for  R and  C were  used  in  Equation  18  and  on  the 
other  the  expressions  for  the  three  "natural"  cooling  processes  were 
also  used  in  the  equation. 
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In  retrospect,  it  can  be  seen  that  the  quantity  called  skin  wettedness 
in  the  new  effective  temperature  model  is  very  similar  to  a heat  stress 
index.  Instead  of  dividing  the  maximum  possible  evaporative  cooling  rate 
into  a required  evaporative  cooling  rate,  the  maximum  value  is  divided 
into  a sweat  cooling  rate  that  is  predicted  from  known  physiological 
responses.  In  line  with  this  thinking,  the  results  of  Gagge*s  model  are 
replotted  in  Figure  22  and  interpreted  as  another  heat  stress  index. 

4.7  Thermal  Comfort  When  Equilibrium  is  Maintained  by  Sweating 

Before  publishing  the  new  effective  temperature  model,  Gagge  and 
his  associates  published  still  another  model  for  indicating  thermal  com- 
fort (18)  that  is  very  similar  to  the  heat  stress  index.  The  cooling 
due  to  evaporation  of  sweat  production  necessary  to  maintain  body  thermal 
equilibrium  is  determined  by  subtracting  the  body  heat  loss  by  radiation, 
convection,  latent  loss  from  the  lung  passages  and  the  loss  due  to  evapora- 
tion of  moisture  diffusing  through  the  skin  from  the  metabolic  heat  pro- 
duction. (See  Appendix  A for  the  specific  relations.)  The  resulting 
required  evaporative  cooling  is  divided  by  the  maximum  possible  and  then 
multiplied  by  100  to  give  a factor  "percentage  of  wettedness". 
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The  main  difference  between  this  model  and  the  two  previous  ones 
mentioned  is  that  the  skin  temperature  is  not  held  constant  as  in  the 
case  of  heat  stress  index  or  predicted  from  known  physiological  responses 
as  with  the  new  effective  temperature  model,  but  rather  assumed  to  be  a 
polynomial  function  of  the  environment  temperature  (see  Appendix  A). 

The  function  was  determined  by  fitting  an  equation  (by  regression 
analysis)  to  experimental  data  that  was  determined  at  the  John  B,  Pierce 
Foundationo  Computations  were  made  using  the  model  and  the  results  are 
shown  in  Figure  23  for  a resting  person  with  the  standard  clothing 
(0„6  do)  subjected  to  a uniform  environment  with  an  air  velocity  of 
20  fpm.  Since  the  model  is  so  similar  in  structure  to  the  heat  stress 
index  and  the  skin  wettedness  factor  of  the  new  effective  temperature 
model,  all  three  are  compared  in  Figure  24.  Since  all  three  use  the 
same  expressions  for  R,  C,  and  the  discrepancies  can  be  attributed 

to  the  three  different  ways  in  which  the  skin  temperature  is  determined 
and  also  to  the  absence  or  inclusion  of  the  "natural"  cooling  terms 
(see  Appendix  A). 
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4.8  Index  of  Thermal  Stress 


Givoni  (2,  19)  has  done  extensive  work  over  the  past  few  years  on 
measuring  and  evaluating  heat  stress.  Having  compared  actual  response 
with  the  predictions  of  many  models  and  found  much  disagreement,  he 
developed  his  own  "index  of  thermal  stress"  (ITS)  which  according  to 
additional  extensive  experiments  conducted  by  Givoni  correctly  predicts 
the  physiological  strain  imposed  on  resting  and  working  people  by  metabolic 
and  environmental  factors.  The  index  is  structured  exactly  like  the 
heat  stress  index  with  two  important  differences.  An  expression  has 
been  included  for  the  radiant  heat  load  due  to  solar  radiation  so  that 
the  index  can  predict  responses  for  people  working  or  resting  outdoors. 

This  is  not  directly  applicable  to  this  study  but  certainly  makes  a 
useful  addition  to  the  field  of  heat  stress  evaluation.  The  other 
difference  is  a factor  called  cooling  efficiency  of  sweating  (f)  and 
is  included  in  the  equation  as  follows: 


The  concept  is  that  where  other  indices  assume  all  the  sweat  produced 
or  required  is  evaporated  and  provides  cooling  for  the  body,  this  may 
not  be  true.  For  example,  some  of  the  evaporation  could  occur  in  the 
clothing  and  not  result  directly  in  cooling  for  the  body.  Calculations 


(19) 


using  Givoni *s  expressions  (see  Appendix  A)  have  been  made  and  are  shown 


in  Figure  25.  The  data  used  was  for  a resting  person  (M  = 


100  kcalv 
hr 


subjected  to  a uniform  environment  with  an  air  velocity  of  .1  m/sec. 


The  clothing  ensemble  chosen  (affects  the  radiation  and  convection 
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expression)  was  very  similar  to  the  Kansas  State  standard  clothing. 

It  should  be  noted  that  the  upper  limit  for  comfort  (according  to 
Givoni)  is  associated  with  an  evaporative  loss  of  60  grams /hr  and  that 
distinct  thermal  discomfort  is  experienced  progressively  as  the  sweat 
rate  is  elevated  above  100  grams /hr. 

Even  though  Givoni  states  that  the  model  predicts  a sweat  rate  that 
is  in  agreement  with  experimentally  measured  rates,  the  results  shown 
in  Figure  25  appear  questionable.  Sweating  is  seen  necessary  at  approxi 
mately  64  The  index  does  show  small  dependence  on  relative  humidity 

at  the  lower  temperatures  and  somewhat  more  dependence  at  the  higher 
temperature  which  agrees  with  the  new  effective  temperature  model. 
However,  it  should  be  remembered  that  what  is  shown  in  the  figure  is 
sweat  rate  and  not  the  ratio  between  the  sweat  rate  and  some  maximum 
possible.  To  put  it  in  a form  comparable  to  the  heat  stress  index, 
values  of  the  above  mentioned  ratio  were  calculated  and  are  shown  in 
Figure  26  and  compared  with  the  original  ITS  values  in  Figure  27. 

4.9  Predicted  Four  Hour  Sweat  Rate 

The  P4SR  index  was  developed  by  McArdle  and  colleagues  during  World 
War  II  at  the  Royal  Naval  Research  Establishment  (20).  Sweat  rate  was 
measured  for  fit,  young  men  dressed  in  shorts  only  or  overall  and  shorts 
when  subjected  to  many  different  combinations  of  temperature,  humidity, 
and  air  velocities  for  four  hour  periods.  The  results  were  presented 
in  the  form  of  nomographs  and  have  been  adapted  to  a psychrometric  chart 
as  shown  in  Figure  28.  The  particular  conditions  for  this  chart  are  a 
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(mean  radiant  temperature  equal  to  air  temperature),  and  an  air  velocity 
of  10  fpm.  As  was  the  case  with  the  original  heat  stress  index,  sweat 
production  is  shown  to  start  only  after  the  temperature  has  reached 
approximately  80  ®F.  Based  on  the  experiments,  a limit  of  3.0  was  es- 
tablished as  the  upper  tolerable  limit  for  daily  exposure. 


In  1957  Lind  and  Hellon  (21)  proposed  the  use  of  the  Oxford  Index 
or  wet-dry  index  to  assess  the  severity  of  hot  climates.  It  is  simply 
a combination  of  the  wet -bulb  and  dry-bulb  temperatures  (°F)  according 
to 


The  extreme  dependence  on  relative  humidity  is  obvious  and  is  shown  in 
Figure  29*  Due  to  its  simplicity  it  has  been  used  quite  often  in  setting 
limits  or  correlating  data  in  extremely  hot  and  wet  conditions  where 
there  is  no  radiant  stress  and  very  low  air  velocity. 


Another  index  very  similar  to  the  Oxford  Index  is  one  called  the 
temperature-humidity  index  defined  by 


4.10  Wet-Dry  Index 


WD  = 0.85  T , + 0.15 

wb  db 


(20) 


4.11  Temperature-Humidity  Index 


THI  =0.4  (T  , + T , ) + 4.8 
ab  wb 


(21) 
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where  all  values  are  in  ®C.  This  index  like  the  last,  attempts  to  combine 
the  effect  of  temperature  and  humidity  in  environments  where  radiant 
stress  and  air  velocity  are  not  important o Figure  30  shows  that  the 
THI  index  is  much  less  dependent  on  humidity  than  the  W-D  index  was . 

As  an  example  of  its  use,  Bridger  and  Helfand  (22)  correlated  death 
rates  during  the  1966  heat  wave  in  St.  Louis  and  found  that  when  the 
24  hour  average  of  the  THI  was  greater  than  27  °C  (81  °F),  there  was 
an  usually  large  number  of  deaths  due  to  heat  stroke. 


The  authors  have  intended  this  part  to  be  simply  a review  and 
presentation  of  important  heat  stress  and  comfort  indices  that  could  be 
integrated  into  the  air  conditioning  criteria  study.  There  is  no  clear 
cut  choice  where  one  index  has  been  universally  accepted  as  correctly 
predicting  physiological  responses  in  all  situations.  On  the  contrary, 
there  is  considerable  disagreement  on  the  relative  merit  of  the  various 
indices.  However  a general  method  of  attack  can  be  suggested.  Indices 
showing  very  slight  dependence  on  humidity  in  the  probable  comfort  range 
of  75  to  80  ®F  (i.e.,  KSU,  PMV,  new  ET)  could  be  used  for  setting  limits 
on  comfort  conditions.  In  addition,  the  ones  showing  strong  dependence 
on  humidity  at  the  higher  temperatures  (i.e.,  new  ET,  W-D,  HSI)  could 
be  used  for  setting  upper  bounds  for  prevention  of  serious  health  hazards. 
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Other  important  factors  should  be  noted.  Most  of  the  indices  are 
valid  for,  or  at  least  have  been  developed  for,  quasi-steady  situations. 
The  subjects  have  been  exposed  to  (or  the  models  simulate  the  subjects 
being  exposed  to)  steady  state  environments  for  periods  of  time  up  to 
four  hours  and  their  response  is  given  accordingly.  Very  seldom  does 
this  occur  in  the  average  home  or  office  building  and  the  technique  of 
extending  these  indices  to  a more  transient  environment  is  questionable. 
ASHRAE,  in  their  comfort  standard,  recommend  limiting  the  rate  of  change 
of  dry  bulb  temperature  in  a space  to  4 °F/hr  and  the  rate  of  change  of 
relative  humidity  to  20%/hour.  These  limits  have  been  verified  somewhat 
by  an  experimental  study  at  Kansas  State  University  reported  by  Sprague 
and  McNall  (23).  College  students  wearing  the  standard  clothing  were 
subjected  to  an  environment  where  the  dry  bulb  temperature  and  relative 
humidity  were  oscillated  one  at  a time  about  the  comfort  value  while  the 
other  parameter  was  maintained  at  the  comfort  value.  The  conclusion  of 
the  study  was  that  the  comfort  standard  values  are  somewhat  conservative 
but  no  change  was  recommended. 
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Another  factor  to  be  considered  quite  important  is  the  effect  of 
nonuniform  environments.  There  is  a capability  built  into  many  of  the 
indices  to  handle  a situation  where  the  wall  temperatures  (all  considered 
the  same)  are  different  than  the  ambient  air  temperature.  However,  one 
parameter  such  as  the  mean  radiant  temperature  cannot  possibly  describe 
actual  environments  where  perhaps  every  surface  is  different  in  tempera- 
ture and  may  not  even  be  uniform  on  the  individual  surfaces.  In  addition, 
very  little  or  nothing  is  known  about  the  physiological  response  to  this 
kind  of  environment.  All  of  the  considerations  mentioned  above  simply 
point  to  the  tremendous  amount  of  work  to  be  done  in  this  field. 

5.  Feasibility  Study  for  the  Determination  of 
an  Air  Conditioning  Criteria 

In  order  to  determine  whether  the  approach  mentioned  in  Part  1 
is  feasible  for  determining  a meaningful  air  conditioning  criteria,  an 
analysis  depicted  in  Figure  31  has  been  carried  out.  An  apartment  to  be 
built  in  Macon,  Georgia  and  Jersey  City,  New  Jersey  as  part  of  HUD*s 
^'Operation  Breakthrough**  (Camci)  has  been  simulated  in  the  National 
Bureau  of  Standards*  computer  program  (NBSLD)  described  in  Part  2. 
Required  as  input  was  all  pertinent  building  data  for  the  apartment  as 
well  as  hour  by  hour  values  of  important  weather  parameters.  Hourly 
values  of  outdoor  dry-bulb  temperature,  wet-bulb  temperature,  wind 
velocity,  and  cloud  cover  were  taken  from  the  United  States  Weather 
Bureau  tapes  (available  from  the  National  Climate  Center  in  Ashville, 

North  Carolina)  for  the  months  of  June,  July,  August  and  September, 

1949  through  1958.  The  output  of  the  program  was  hourly  values  of 
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indoor  dry-bulb  and  wet-bulb  temperature  for  the  four  month,  ten  year 
period.  Algorithms  for  pertinent  physiological  indices  were  written 
along  with  detailed  statistical  subroutines  so  that  daily  profiles  and 
certain  histograms  showing  values  of  these  indices  could  be  obtained. 
The  building  data  that  was  used  is  as  follows: 

1.  The  apartment  building  was  assumed  to  face  west  and  the 
individual  apartment  chosen  was  one  at  approximately 
mid-height  of  the  130  feet  high  structure.  No  building 
shadow  was  used  and  the  adjacent  apartments  and  hallway 
were  assumed  to  always  be  at  the  same  temperature  as 
the  apartment . 

2.  The  exposed  wall  is  to  be  constructed  of  7 inches  of 

concrete  and  2 inches  of  insulation.  It  contains  a 
2 

100  ft  glass  window  (shading  coefficient  = 0.22 
corresponding  to  a double  pane  window  with  a white 

2 

opaque  roller  shade).  The  total  wall  area  is  212  ft 
2 

leaving  112  ft  of  solid  structure. 

3.  The  floor/ceiling  unit  is  composed  of  5 1/2  inches  of 

concrete  and  1/8  inch  cork  tile  with  a total  area  of 
2 

550  ft  . The  floor  to  ceiling  height  is  8.5  ft. 

4.  The  partitioned  walls  contain  6 inches  of  concrete. 

5.  The  maximum  number  of  occupants  is  2 , the  maximum 

2 

equipment  load  is  1/2  watts /ft  of  floor  area,  and 

2 

the  maximum  lighting  load  is  3 watts/ft  of  floor  area. 
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6,  The  infiltration  or  natural  ventilation  rate  was 
simulated  as  1 air  change/hr  between  6 a.m,  and  7 
p.m.  and  6 air  changes/hour  from  7 p.m,  to  6 a.m. 

Typical  results  for  the  four  months  of  one  summer  in  Jersey  City 
are  shown  in  Figures  32  through  67.  The  data  is  for  the  summer  of  1954 
which  was  perhaps  the  coolest  of  the  10  year  period.  The  wet -bulb 
temperature  was  above  67  °F  for  only  583  hours  and  the  dry-bulb  tem- 
perature was  above  80  °F  for  only  493  hours  of  the  four  month  summer. 
Shown  in  Figure  32  is  a plot  of  the  solar  energy  falling  on  the  apart- 
ment (after  modification  by  actual  cloud  cover)  and  a plot  of  outdoor 
dry  and  wet-bulb  temperature  for  the  month  of  June.  Figure  33  shows 
the  output  of  the  thermal  simulation  of  this  apartment;  in  other  words, 
a monthly  profile  of  indoor  dry  and  wet  bulb  temperature.  Since  the 
ASHRAE  comfort  standard  gives  limits  on  the  rate  of  change  of  indoor 
temperature  and  relative  humidity,  these  values  were  computed  and  are 
shown  in  Figure  34.  As  can  be  seen,  the  suggested  values  of  4 °F/hr 
and  20%/hr  were  only  exceeded  once  or  twice  during  the  month  and  this 
was  also  found  true  for  every  other  summer  month  of  the  ten  year  period 
in  Jersey  City  and  Macon,  regardless  of  how  hot  the  weather  was.  Figures 
35  through  40  show  monthly  profiles  of  several  of  the  indices  discussed 
in  Part  4.  Comments  that  can  be  made  concerning  these  profiles  are 
as  follows: 
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1. 


The  new  effective  temperature  was  higher  than  the  upper 
limit  of  the  comfort  zone  (81)  8,7  percent  of  the  time 
and  higher  than  the  upper  limit  of  the  broad  comfort 
zone  (85)  only  1 percent  of  the  time  during  the  month. 

2.  The  predicted  mean  vote  for  a person  in  standard 
clothing  (0.6  do)  exceeded  1 (slightly  warm)  4.2 
percent  of  the  time  but  was  above  1 only  1.1  percent 
of  the  time  when  the  clothing  was  reduced  (0.25  do). 

3.  There  was  practically  little  or  no  heat  stress 
regardless  of  the  particular  index  used  (old  or  new). 

4.  The  Kansas  State  Index  was  above  5 (slightly  warm) 
only  1.0  percent  of  the  time. 

5.  More  than  10  percent  of  the  people  were  dissatisfied 
only  four  times  during  the  month  with  light  clothing 
but  that  increased  to  7 or  8 times  with  standard 
clothing  assumed. 

The  same  kinds  of  observations  could  be  made  about  the  three  remaining 
months  from  Figures  41  through  67. 
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In  addition  to  the  summer  profiles  of  the  various  parameters, 
certain  statistical  data  was  gathered  and  examples  are  shown  in  Table 
5 and  Figures  68  through  70,  The  table  shows  for  the  10  summers  in 
Jersey  City,  the  number  of  total  hours  in  each  summer  that  the  outdoor 
dry  bulb  temperature  exceeded  80  °F  and  the  outdoor  wet -bulb  temperature 
exceeded  67  °F.  In  addition,  the  percent  of  the  time  in  each  month 
that  selected  indices  exceeded  prescribed  levels  was  calculated.  Since 
not  only  the  total  length  of  time  but  also  duration  of  each  hot  period 
will  affect  the  general  satisfaction  or  dissatisfaction  with  the  space 
conditions,  the  length  of  time  that  each  index  stayed  above  the  pre- 
scribed level  for  each  hot  period  was  computed.  In  Figure  68  for 
example,  the  new  effective  temperature  was  above  85  for  a total  of  14 
percent  of  the  time  during  August  of  1953,  The  specific  way  it  occurred 
was : once  for  2 hours , twice  for  between  12  and  18  hours  and  once  for 

more  than  48  hours.  Contrast  this  with  July  of  1958  (Figure  69)  where 
the  total  amount  of  time  was  approximately  the  same  but  the  manner  in 
which  it  happened  was  entirely  different:  4 times  for  only  2 hours  at 

a time,  5 times  for  3 hours,  once  for  6 hours,  twice  for  7 hours,  twice 
for  8 hours,  once  for  10  hours,  and  finally  once  for  between  12  and  18 
hours.  Most  probably,  persons  would  not  express  the  same  degree  of 
dissatisfaction  during  these  two  months. 


54 


It  is  felt  that  given  the  kind  of  information  just  shown  for  a 
specific  structure  and  location  the  decision  to  air  condition  or  not 
could  certainly  be  made.  However  it  would  be  neither  practical  nor 
economical  to  undertake  such  an  analysis  for  each  individual  case. 

It  may  not  even  be  feasible  to  carry  out  the  analysis  for  all  possible 
combinations  of  building  groups  and  climatic  zones  across  the  United 
States.  Therefore  a short-cut  method  was  sought  by  examining  the 
results  of  this  thermal  simulation  very  closely. 

In  Figure  71  through  76  it  can  be  seen  that  correlations  were 
attempted  between  time  (in  4 months)  above  a prescribed  index  value 
inside  and  time  the  outdoor  dry-bulb  or  wet-bulb  temperature  exceeded 
80  °F  or  67  °F  respectively.  Even  if  successful,  the  result  would 
naturally  have  only  been  valid  for  the  specific  apartment  in  the  specific 
location.  However  it  was  felt  that  this  would  be  one  way  to  describe 
quantitatively  the  "thermal  performance"  of  the  unit.  As  the  figures 
indicate,  the  number  of  hours  in  several  consecutive  months  that  the 
dry-bulb  temperature  is  above  80  ®F  or  the  wet-bulb  temperature  is  above 
67  °F  does  not  correlate  well  with  what  occurs  inside  the  apartment 
during  this  same  period  of  time.  These  parameters  of  course  are  pre- 
cisely the  ones  used  in  the  present  HUD  criteria. 
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Seeking  a correlation  covering  a shorter  period  of  time,  Figures 
77  and  78  show  the  daily  maximum  values  of  the  KSU  index  inside  as  a 
function  of  maximum  dry-bulb  temperature  and  wet -bulb  temperature  outside 
for  August  of  1954  in  Jersey  City.  The  results  are  much  better  than  the 
previous  ones  but  still  too  scattered  to  be  of  much  value.  However,  when 
the  same  dependent  variable  was  plotted  as  a function  of  average  (over 
a 24  hour  period)  outdoor  dry-bulb  temperature  the  result  was  excellent 
as  shown  in  Figure  79.  The  correlation  coefficient  (r)  for  this  particu- 
lar set  of  data  was  0.97185  (maximum  possible  value  = 1.0).  The  same 
degree  of  success  was  attained  with  other  indices  such  as  the  new  ef- 
fective temperature  shown  in  Figure  80.  The  daily  average  temperature 
used  was  not  the  average  of  the  twenty-four  hourly  values  but  rather 
the  average  of  the  daily  maximum  and  daily  minimum.  This  particular 

average  was  chosen  due  to  the  fact  that  many  United  States  Weather  Bureau 

* / 

Stations  report  precisely  this  parameter"  . The  discrepancy  is  usually 
quite  small  as  the  curves  of  Figure  81  indicate. 


~ Instructions  for  Climatological  Observers  Circular  B,  U.  S.  Weather 
Bureau  Publication,  January,  1962. 
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The  straight  lines  of  Figure  79  and  80  were  determined  by  regression 
analysis  (method  of  least  squares)  in  order  to  fit  the  data  as  close  as 
possible.  In  fact,  similar  regression  lines  were  determined  for  June, 

July  and  August  of  all  ten  years  in  both  Jersey  City  and  Macon.  It  was 
found  that  each  of  the  30  lines  for  Jersey  City  fell  within  the  corridor 
of  Figure  82  and  similarly  in  Macon  (Figure  83).  As  would  be  expected 
the  changing  of  any  major  building  parameter  would  cause  the  apartment 
to  "perform"  differently.  The  effect  of  only  reducing  the  shading  over 
the  large  window  can  be  seen  in  Figure  84  and  similarly  in  Figure  85  for 
keeping  everything  the  same  except  the  facing  direction  of  the  exterior  wall« 

The  importance  of  such  correlation  results  lie  in  the  fact  that  if 
proven  valid  by  additional  analytical  and/or  experimental  studies  a 
particular  line  could  be  determined  for  a given  housing  unit-weather 
data  group  by  a small  number  of  calculations  compared  to  a 4 month  10 
year  hour  by  hour  computation.  The  decision  to  air  condition  would  then 
have  to  be  made  not  on  the  total  time  or  percent  of  time  over  an  extended 
period  that  an  index  value  was  exceeded  but  rather  whether  the  daily 
average  outdoor  temperature  for  some  given  design  conditions  caused  an 
indoor  environment  that  was  unacceptable  when  compared  with  prescribed 
levels  of  the  PIHI. 
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To  compare  an  approach  such  as  this  with  the  ten  year  computation 
already  presented,  additional  calculations  were  carried  out.  The  1% 
design  values  of  outdoor  dry-bulb  and  wet-bulb  temperature  presented 
in  the  ASHRAE  Handbook  of  Fundamentals  (24)  for  several  selected  cities 
(see  Table  6)  were  used  in  conjunction  with  the  dimensionless  daily  tem- 
perature cycle  of  Figure  86  (determined  from  surveying  actual  weather 
tapes)  to  produce  a design  day  cycle  for  these  cities.  The  Camci  apart- 
ment was  then  subjected  to  this  design  day  cycle  (in  the  computer  program) 
for  seven  straight  days  using  the  actual  solar  radiation  in  these  locali- 
ties and  assuming  zero  cloud  cover.  Shown  in  Figure  87  are  the  results 
of  the  calculations  where  the  maximum  value  of  the  KSU  index  in  the  apart- 
ment during  the  second  and  fifth  day  are  plotted  for  each  city.  There 
was  no  noticeable  change  after  the  fifth  day.  As  can  be  seen,  the 
agreement  between  the  data  points  and  previously  determined  corridors 
is  quite  good. 


58 


6 .  References 


1.  J.  F.  van  Straaten,  "Thermal  Performance  of  Buildings",  Elsevier 
Publishing  Company,  1967. 

2.  B.  Givoni,  "Man,  Climate,  and  Architecture",  Elsevier  Publishing 
Company,  1969. 

3.  R.  G.  Kevins  and  P.  E.  McNall,  Jr.,  "Proposed  Classification  of 
Human  Response  to  the  Environment",  ASHRAE  Journal,  pp.  44-47, 

June  1971. 

4.  Procedure  for  Determining  Heating  and  Cooling  Loads  for  Computerized 
Energy  Calculations;  Algorithms  for  Building  Heat  Transfer  Subroutines, 
ASHRAE,  1971. 

5.  T.  Kusuda,  "Thermal  Response  Factors  for  Multilayer  Structures  of 
Various  Heat  Conduction  Systems",  ASHRAE  Transactions,  Vol.  75, 

Part  1,  1969. 

6.  F.  C.  Houghten  and  C.  P.  Yaglou,  "Determination  of  the  Comfort  Zone", 
ASHVE  Transactions,  Vol.  29,  p.  361,  1923. 

7.  T.  Bedford,  "Environmental  Warmth  and  Its  Measurement",  MRC  War 
Memorandum  No.  17,  HMSO,  London,  1946. 

8.  N.  Koch,  B.  H.  Jennings  and  C.  M.  Humphreys,  "Sensation  Responses 

to  Temperature  and  Humidity  Under  Still  Air  Conditions  in  the  Comfort 
Range",  ASHRAE  Transactions,  Vol.  66,  p.  264,  1960. 

9.  H.  Missenard,  "Equivalences  thermiques  des  ambiences;  equivalences 
de  passage;  equivalences  de  sejour",  Chaleur  et  Industrie,  July- 
August,  1948. 


59 


10 


. R.  G,  Nevins,  F.  H.  Rohles,  W.  E,  Springer,  and  A.  M.  Feyerherm, 

*'A  Temperature  Humidity  Chart  for  Thermal  Comfort  of  Seated  Per- 
sons", ASHRAE  Transactions,  Vol.  72,  Part  1,  1966, 

11.  R.  G.  Nevins  and  F.  H.  Rohles,  "The  Nature  of  Thermal  Comfort  for 
Sedentary  Man",  ASHRAE  Transactions,  Vol.  77,  Part  1,  1971. 

12.  0.  Seppanen,  P.  E.  McNall,  D.  M.  Munson,  and  C.  H.  Sprague,  "Thermal 
Insulating  Values  for  Typical  Clothing  Ensembles",  to  be  published. 

13.  F.  H.  Rohles,  Jr.,  "The  Modal  Comfort  Envelope  - A New  Approach 
Toward  Defining  the  Thermal  Environment  in  which  Sedentary  Man  is 
Comfortable",  ASHRAE  Transactions,  Vol.  76,  Part  1,  1970. 

14.  P.  0.  Fanger,  "Thermal  Comfort",  Danish  Technical  Press,  Copenhagen, 
1970. 

15.  A.  P.  Gagge,  J.  A.  J.  Stolwijk,  and  Y.  Nishi,  "An  Effective  Tem- 
perature Scale  Based  on  a Simple  Model  of  Human  Physiological 
Regulatory  Response",  ASHRAE  Transactions,  Vol.  77,  Part  1,  1971. 

16.  H.  S.  Belding  and  T.  F.  Hatch,  "Index  for  Evaluating  Heat  Stress 
in  Terms  of  Resulting  Physiological  Strain",  Heating,  Piping,  Air 
Conditioning,  Vol.  27,  p.  129,  1955. 

17.  J.  S.  McKarns,  and  R.  S.  Brief,  "Nomographs  Give  Refined  Estimates 
of  Heat  Stress  Index",  Heating,  Piping,  and  Air  Conditioning,  Vol. 
38,  p.  113,  1966. 

18.  A.  P.  Gagge,  J.  A.  J.  Stolwijk,  and  Y.  Nishi,  "The  Prediction  of 
Thermal  Comfort  When  Equilibrium  is  Maintained  by  Sweating",  ASHRAE 
Transactions,  Vol.  75,  Part  II,  p.  108,  1969. 


60 


19.  B.  Givoni,  "Estimation  of  the  Effect  of  Climate  On  Man:  Develop- 

ment of  a New  Thermal  Index",  Research  Report  to  UNESCO,  Building 
Research  Station,  Technion,  Haifa,  Israel,  1963. 

20.  Be  McArdle,  W,  Durham,  H.  E.  Rolling,  W.  S.  S.  Ladell,  J.  W.  Scott, 

M.  L.  Thomson,  and  J.  S.  Weiner,  "The  Prediction  of  the  Physiological 
Effects  of  Warm  and  Hot  Environments",  MRC  R.N.P.,  Vol.  47,  p.  391, 
1947. 

21.  A.  R.  Lind,  and  R.  F.  Hellon,  "Assessment  of  Physiological  Severity 
of  Hot  Climates",  Journal  of  Applied  Physiology,  Vol.  II,  p.  35, 

1957  c 

22.  C.  A.  Bridger  and  L.  A.  Helfand,  "Mortality  from  Heat  During  July 
1966  in  Illinois",  Int.  Jour.  Biomet.,  Vol.  12,  p.  51,  1968. 

23.  C.  H.  Sprague,  and  P.  E.  McNall,  Jr.,  "The  Effects  of  Fluctuating 
Temperature  and  Relative  Humidity  on  the  Thermal  Sensation  (Thermal 
Comfort)  of  Sedentary  Subjects",  ASHRAE  Transactions,  Vol.  76, 

Part  1,  1970. 

24.  ASHRAE  Handbook  of  Fundamentals,  ASHRAE,  1967. 

25.  T.  H.  Benzinger,  "Clinical  Temperature,  New  Physiological  Basis", 
Journal  of  the  American  Medical  Association,  Vol.  209,  pp.  1200- 
1206,  1969. 

26.  T.  H.  Benzinger,  "Peripheral  Cold  Reception  and  Central  Warm 
Reception,  Sensory  Mechanisms  of  Behavioral  and  Autonomic  Thermo - 
stasis".  Physiological  and  Behavioral  Temperature  Regulation, 

Chapter  56,  Charles  C.  Thomas,  Springfield,  Illinois,  1970. 


61 


27.  J.  D.  Hardy  and  J.  A.  J.  Stolwijk,  "Partitional  Calorimetry  Studies 
of  Man  During  Exposure  to  Thermal  Transients",  Journal  of  Applied 
Physiology,  Vol.  21,  pp.  1799-1806,  1966. 

28.  J.  A.  Jo  Stolwijk  and  J.  D.  Hardy,  "Temperature  Regulation  in  Man  - 
A Theoretical  Study",  Pfltigers  Archiv.,  Vol.  291,  pp.  129-162,  1966. 

29.  T.  H.  Benzinger  and  C.  Kitzinger,  "Gradient  Layer  Calorimetry  and 
Human  Calorimetry",  Temperature  - Its  Measurement  and  Control  in 
Science  and  Industry,  Vol.  3,  Part  3,  Reinhold  Publishing  Corp., 

N.  Y.,  N.  Y.,  1963. 


62 


Appendix  A 


Component  Expressions  for  Various  Comfort  Models 
and  Ranges  of  Application 

In  Part  4,  the  methods  of  computation  for  several  of  the  analyti- 
cally determined  comfort  indices  were  discussed.  Since  most  of  them  use 
expressions  for  heat  transfer  rate  between  the  subject  and  environment 
by  convection,  radiation,  etc.,  these  expressions  are  given  here  for 
purposes  of  comparison. 

Rate  of  Heat  Loss  by  Convection  (C) 

PMV:  C = A^„  f , h (T  1 - T ) kcal/hr 
dU  c 1 c cl  a 

h = 2.05  (T  - - T for  2.05  (T  , - T > 10.4 

c c -L  a c X a 

= 10.4  Jv  for  2.05  (T  , - T < 10.4  > 

c X a 

2 

ET  (new):  C=h  (T,-T)F-  watts/m 

c s X a c X 

h = 11.6  v^*^  watts/m^  °C 
c 

HSI  (original):  C = 2 v®'^  (95.  - T ) Btu/hr 

a 

where  v = air  velocity,  fpm 

= ambient  air  temperature,  °F 

HSI:  identical  with  ET  (new) 
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Thermal  Comfort  when  Equilibrium  is  Maintained  by  Sweating:  identical 

with  ET  (new) 

ITS:  C = a (35  - T ) kcal/hr 

a 

Rate  of  Heat  Loss  by  Radiation  (R) 

PMV:  R = 4.8  X 10’®  Ajjjj  [(T^^  + 273)^  - + 273)"^]  kcal/hr 

ET  (new) : R=h  (T,-T)F.  watts/m^ 

r SK  a ci 

HSI  (original):  R = 22  (95  - T ) Btu/hr 

“ mean  radiant  temperature,  °F 

NKl 

HSI:  identical  with  ET  (new) 

Thermal  Comfort  when  Equilibrium  is  Maintained  by  Sweating:  identical 

with  ET  (new) 

ITS:  accounted  for  in  the  convection  term  by  replacing  T^  by 

Rate  of  "Dry”  Heat  Loss  in  Respiratory  Passages  (E  ) 

res 

PMV:  E = .0014  M (34  - T ) kcal/hr 

res  a 
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ET  (new):  zero 


HSI  (original):  zero 

2 

HSI:  E = .0014  M (34  - T ) watts/m 

res  a 

2 

M = metabolic  heat  production,  watts/m 
Thermal  Comfort  when  Equilibrium  is  Maintained  by  Sweating:  zero 


ITS:  zero 


Rate  of  Energy  Loss  due  to  Evaporative  Cooling  in  Respiratory  Passages 


(E  ) 
res 


PMV:  E = .0023  M (44  - P ) kcal/hr 

res  a 


ET  (new):  E = .0023  M (44  - P ) watts/m^ 

res  a 

2 

M = metabolic  heat  production,  watts/m 


HSI  (original):  zero 


HSI:  identical  with  ET  (new) 


Thermal  Comfort  when  Equilibrium  is  Maintained  by  Sweating:  identical 

with  ET  (new) 
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ITS : zero 


Rate  of  Energy  Loss  due  to  Evaporative  Cooling,  Normal  Skin  Diffusion 


PMV 


M 


: “ 0,061  P^)  kcal/hr 

du 


ET  (new):  E,  ...=  0.06  • 2,2  h (P  , - P ) F ^ watts/m*^ 

cLx.tr  c sK  a per 


= convection  coefficient  given  above,  watts/m  °C 


HSI  (original):  zero 


HSI:  identical  with  ET  (new) 


Thermal  Comfort  when  Equilibrium  is  Maintained  by  Sweating: 


= .35  (P  , - P ) F 1 watts/m*^ 
diff  sk  a pci 


ITS:  zero 

Rate  of  Energy  Loss  due  to  Evaporative  Cooling.  Regulatory  Sweat  Pro 

duct ion  (E  , E ) 

sw  rsw 

PMV:  E = .42  (t^^ 50)  kcal/hr 

sw  "DU 

ET  (new):  E = 70  (T  - 36.6)  (T  , - 34.1)  watts/m^ 

rsw  cr  SK 
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for  T ^ 36.6  °C  and  T , ^ 34.1  °C 
cr  sk 

E = 0.0 
rsw 

for  T < 36.6  or  T ,<  34.1  °C 
cr  sk 

HSI  (original):  not  applicable 

HSI:  not  applicable 

Thermal  Comfort  when  Equilibrium  is  Maintained  by  Sweating:  not  applicable 

ITS:  not  applicable 

Maximum  Evaporative  Cooling  Rate  Possible  (E^^^) 

PMV:  not  specified 

2 

ET  (new):  E =2.2  h (P.-P)F  ^ watts/m 

max  c sk  a pci 

2 

h^  = convection  coefficient  given  previously,  watts/m  °C 

4 

HSI  (original):  E = 10  v*  (P  , - P ) Btu/hr 

max  sk  a 

V = air  velocity,  fpm 
HSI:  same  as  ET  (new) 
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Thermal  Comfort  when  Equilibrium  is  Maintained  by  Sweating:  same  as 

ET  (new) 


ITS:  E = P v*^  (42  - P ) kcal/hr 

max  a 


Rate  of  Heat  Conduction  Through  Clothing  (K) 

This  particular  quantity  is  only  applicable  to  PMV  and  is  given  by: 


K = A_^.  ^sk  ^cl  kcal/hr 

.18  I T 
cl 


Relation  Between  PMV  and  Body  Heat  Load  (L) 


PMV 


= [.352  + .032]  L 


Relation  for  Sweating  Efficiency 


1 _ 0.6  (E  /E  - 0.12) 

“ “ e req  max 


Relation  for  PPD 

A second  order  equation  was  fit  to  the  tabular  data  given  by 
Fanger  for  the  Predicted  Percentage  of  Dissatisfied  (due  to  heat  only) 
as  a function  of  PMV 


PPD  = 1.63  + 10.98  (PMV)  + 13.41  (PMV)' 
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Relation  for  Skin  Temperature  as  a Function  of  Ambient  Air  Temperature 


The  following  equation  was  used  in  Gagge*s  "Thermal  Comfort  When 
Equilibrium  is  Maintained  by  Sweating"  model 


T , = 25.49  + .249  T °F 
sk  a 


Unless  it  was  otherwise  specified,  the  following  nomenclature 
applies  to  the  relations: 

2 

= body  surface  area,  m 

= clothing  coefficient  (approximately  11.6  for  0.6  do 
standard  clothing) 

f T = ratio  of  the  surface  area  of  the  clothed  body  to  the 
cl 

nude  body  (=  1.1  for  0.6  do  standard  clothing) 

= a factor  that  measures  the  efficiency  for  the  passage 

of  dry  heat  from  the  skin  surface  through  the  clothing 

to  the  environment  (=  .58  for  0.6  do  clothing) 

f _-  = the  ratio  of  the  effective  radiation  area  of  the 
eff 

clothed  body  to  the  surface  area  of  the  clothed 
body  (0.65  for  seated  position  and  0.75  for  standing) 
Fp^^  = permeation  efficiency  factor  for  water  vapor  evaporated 
from  the  skin  surface  through  clothing  to  the  ambient 
air  (=  0.82  for  0.6  do  clothing) 

2 

h = convection  coefficient,  kcal/hr  m °C 
c 
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M = 
P = 


T 


cr 


2 

radiation  coefficient,  watts/m  ®C,  **=  5.23  in  the  ET  model" 
conduction  resistance  value  for  clothing,  do 
metabolic  heat  production,  kcal/hr 

clothing  coefficient  (approximately  13.0  for  standard 
clothing) 

vapor  pressure  in  ambient  air,  mm  Hg 

saturated  vapor  pressure  corresponding  to  skin  temperature, 
mm  Hg 

ambient  air  temperature,  °C 
clothing  surface  temperature,  ®C 
core  or  deep  body  temperature,  ®C 
mean  radiant  temperature,  ®C 
skin  surface  temperature,  ®C 
air  velocity,  m/sec 


In  addition  to  the  relations.  Table  7 has  been  included  to  give 
the  range  of  application  of  the  various  indices  as  given  on  psychrometric 
charts  in  this  report. 
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Appendix  B 


Discussion  of  Regulatory  Sweating  Models 

In  Part  4,  it  was  noted  that  the  component  expression  for  regula- 
tory sweating  in  the  new  effective  temperature  model  contributed  signifi 
cantly  to  the  response  of  the  model.  As  given  in  Appendix  A,  the  com- 
ponent expression  is 

E = 70  (T  - 36.6)  (T  , - 34.1)  watts/m^ 
rsw  cr  sk 

for  T ^36.6  °C  and 
cr 

T , ^34.1  °C 
sk 

E = 0 if  T < 36.6  °C  or 
rsw  cr 

T , < 34.1  °C 
sk 

The  interpretation  of  this  expression  is  that  sweat  production  is  caused 
by  signals  from  both  the  skin  and  the  deep  body.  If  the  temperature  in 
either  one  of  the  regions  rises  above  the  normal  value,  then  the  body 
secrets  sweat  for  the  purpose  of  evaporative  cooling.  This  sweating 
model  is  in  definite  constrast  to  the  results  of  experiments  conducted 
by  Dr.  T.  H.  Benzinger  of  the  Naval  Medical  Research  Institute  (25,  26). 
Dr.  Benzinger  concludes  as  a result  of  his  studies,  that  the  signal  for 
sweat  secretion  is  a sole  function  of  the  temperature  in  the  anterior 
portion  of  the  hypothalamus  region  of  the  brain.  When  this  temperature 
rises  above  a sharply  determined  set  point,  warm  sensitive  neurons  begin 
to  fire  and  excite  sweating  such  that  the  rate  of  sweat  production  is 
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directly  proportional  to  deviation  beyond  the  set  point. 

The  skin  temperature  does  play  a role  in  the  evaporative  cooling 
process;  however,  in  contrast  to  the  interpretation  of  Hardy  and  Stolwijk 
(27,  28)  (which  has  been  adapted  for  the  new  effective  temperature  model), 
Benzinger  concludes  that  the  effect  is  one  of  inhibition.  Whenever  the 
skin  temperature  drops  below  approximately  33  ®C  the  rate  of  sweat  pro- 
duction begins  to  decrease  below  its  normal  value  (given  a hypothalamus 
temperature  above  the  set  point)  and  the  amount  of  decrease  is  proportional 
to  the  deviation  of  skin  temperature  below  33  ®C.  A comparison  between 
these  two  distinctly  different  interpretations  is  shown  in  Figure  88. 

The  rate  of  evaporative  cooling  is  plotted  against  the  deviation  of 
skin  temperature  from  33  °Co  As  can  be  seen,  Gagge*s  model  simulates 
zero  cooling  by  sweat  evaporation  as  long  as  the  skin  temperature  is 

below  34.1  ®C  and  is  directly  proportional  to  (T  , - 34.1)  • (T  - 36.6) 

s K cr 

above  this  point.  Benzinger *s  data  indicates  a constant  rate  of  sweat 
production  whenever  the  skin  temperature  is  above  33  ®C  (the  value  of 
the  rate  dependent  on  (T^^  “ 36.6))  and  a progressively  decreasing  rate 
as  the  skin  temperature  drops  below  33  ®C. 
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The  importance  of  the  evaporative  cooling  in  the  model  can  be  seen 
in  Figure  89  where  new  effective  temperature  lines  are  plotted  that  have 
been  calculated  using  Benzinger*s  results  for  sweat  production.  Since 
his  data  generally  gives  larger  values  of  evaporative  cooling  for 
specified  core  and  skin  temperatures,  it  would  be  expected  that  this 
modified  scale  would  show  more  dependence  on  humidity  than  the  standard 
new  effective  temperature  lines  and  this  can  be  seen  in  Figure  90.  It 
is  not  recommended  that  the  new  data  be  used  in  the  simulation  of  the 
body  response  since  the  experiments  used  to  generate  them  were  very 
special  in  nature.  Men  were  enclosed  in  a gradient  calorimeter  (see 
reference  29)  such  that  all  sweat  evaporated  could  be  instantly  carried 
away.  This  may  not  be  true  in  ordinary  environments,  particularly  where 
the  humidity  of  the  ambient  air  is  high.  However,  it  should  be  obvious 
that  much  more  study  and  refinement  is  needed  before  an  accurate  simula- 
tion of  human  response  to  an  environment  can  be  accomplished. 

The  interpretation  that  lines  of  equivalent  thermal  comfort  are 
identical  with  lines  of  constant  sweat  ratio  could  be  adjusted  in  light 
of  Dr.  Benzinger*s  findings.  The  feeling  of  thermal  comfort  is  no  doubt 
linked  to  the  necessity  for  sweat  secretion  when  subjected  to  warm 
environments.  However,  as  a result  of  his  studies  this  is  in  turn  linked 
directly  to  the  hypothalamus  temperature  or  deep  body  temperature  since 
this  is  the  sole  source  of  the  sweating  signal.  The  reasoning  is  sub- 
stantiated by  an  additional  experiment  conducted  by  him.  A subject  was 
cooled  in  a tank  of  cold  water  until  the  hypothalamus  temperature  was 
well  below  the  individual’s  set  point.  He  was  then  put  into  a very  hot 
bath  where  the  skin  temperature  and  hypothalamus  temperature  were  monitored 
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with  time.  Even  though  the  skin  temperature  went  to  a relatively  high 
value  almost  immediately,  the  subject  did  not  report  discomfort  until 
the  hypothalamus  temperature  began  to  rise  above  his  set  point. 

The  new  effective  temperature  model  as  presently  structured  was 
used  to  generate  lines  of  constant  core  temperature  and  these  are  shown 
in  Figure  91.  The  numerical  values  assigned  to  the  lines  were  again 
determined  by  the  dry-bulb  temperature  along  the  507o  relative  humidity 
curve • As  would  be  expected,  the  lines  don't  show  quite  as  much  de- 
pendence on  humidity  as  before  and  this  can  be  seen  in  Figure  92 « How- 
ever, the  agreement  with  the  studies  at  Kansas  State  University  are 
quite  good  as  shown  in  Figure  93.  In  fact , the  85  ET  line  which  was 
recommended  as  the  upper  limit  of  the  broad  comfort  zone  is  identical 
to  the  KSU  = 5 (slight ly  warm)  line. 
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area  RATiOf  EXT.WALL/FlOOR  = .B59 
area  KATlOt  a/INDOW/EXT.wALL  = .157 


Table  4 


Evaluation  of  Index  of  Heat  Stress 


Index  of 

Heat  Physiological  and  Hygienic  Implications  of  8-Hr  Ex- 

Stress  posures  to  Various  Heat  Stresses 


0 


No  thermal  strain 


+10 


40 

50 

60 


Mild  to  moderate  heat  strain.  Where  a job  involves 
higher  intellectual  functions,  dexterity,  or  alert- 
ness, subtle  to  substantial  decrements  in  perform- 
ance may  be  expected.  In  performance  of  heavy  physi- 
cal work,  little  decrement  expected  unless  ability  of 
individuals  to  perform  such  work  under  no  thermal 
stress  is  marginal. 

Severe  heat  strain,  involving  a threat  to  health 
unless  men  are  physically  fit.  Break-in  period 
required  for  men  not  previously  acclimatized.  Some 
decrement  in  performance  of  physical  work  is  to  be 
expected.  Medical  selection  of  personnel  desirable 
because  these  conditions  are  unsuitable  for  those 
with  cardiovascular  or  respiratory  impairment  or  with 
chronic  dermatitis.  These  working  conditions  are 
also  unsuitable  for  activities  requiring  sustained 
mental  effort. 


70 

80 

90 


100 


Very  severe  heat  strain.  Only  a small  percentage  of 
the  population  may  be  expected  to  qualify  for  this 
work.  Personnel  should  be  selected  (a)  by  medical 
examination  and  (b)  by  trial  on  the  job  (after 
acclimatization).  Special  measures  are  needed  to 
assure  adequate  water  and  salt  intake.  Amelioration 
of  working  conditions  by  any  feasible  means  is  highly 
desirable,  and  may  be  expected  to  decrease  the  health 
hazard  while  increasing  efficiency  on  the  job.  Slight 
"indisposition"  which  in  most  jobs  would  be  insufficient 
to  affect  performance  may  render  workers  unfit  for  this 
exposure o 

The  maximum  strain  tolerated  daily  by  fit,  acclimatized 
young  men. 


Month  Total  No.  of  Hours  Total  No.  of  Hours  Percent  of  Percent  of  Percent  of 
and  in  4 Months  in  4 Months  Time  New  Time  New  Time  KSU 
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Table  6 ASHRAE  1%  Design  Weather 
Conditions  for  Selected  Cities 


Location 

^db> 

Range , °F 

Range 
db  2 

Portland,  Maine 

88.0 

22.0 

77.0 

Concord,  N,  H. 

91.0 

26.0 

78.0 

Har  t f or d , Conn  o 

90.0 

22.0 

79.0 

West  Chester,  Pa» 

92.0 

20.0 

82.0 

New  York,  N.  Y. 
(Kennedy ) 

91.0 

16.0 

83.0 

Jersey  City,  N.  J. 

94.0 

20.0 

84.0 

Washington,  D,  Co 

94.0 

18.0 

85.0 

New  York,  N,  Y. 
(Central  Park) 

94.0 

17.0 

85.5 

Table  7 Range  of  Application  of  Various  Physiological  Indices  as  Calculated  in  This  Report 
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